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ABSTRACT 

Intense and complex instrumental backgrounds, against which the much smaller sig- 
nals from celestial sources have to be discerned, are a notorious problem for low and 
intermediate energy 7 -ray astronomy (~ 50 keV - 10 MeV). Therefore a detailed qualita- 
tive and quantitative understanding of instrumental line and continuum backgrounds is 
crucial for most stages of 7 -ray astronomy missions, ranging from the design and devel- 
opment of new instrumentation through performance prediction to data reduction. We 
have developed MGGPOD, a user-friendly suite of Monte Carlo codes built around the 
widely used GEANT (Version 3.21) package, to simulate a b initio the physical processes 
relevant for the production of instrumental backgrounds. These include the build-up 
and delayed decay of radioactive isotopes as well as the prompt de-excitation of excited 
nuclei, both of which give rise to a plethora of instrumental 7 -ray background lines in 
addition to continuum backgrounds. The MGGPOD package and documentation are 
publicly available for download. 

We demonstrate the capabilities of the MGGPOD suite by modeling high resolution 
7 -ray spectra recorded by the Transient Gamma-Ray Spectrometer (TGRS) on board 
Wind during 1995. The TGRS is a Ge spectrometer operating in the 40 keV to 8 MeV 
range. Due to its fine energy resolution, these spectra reveal the complex instrumental 
background in formidable detail, particularly the many prompt and delayed 7 -ray lines. 
We evaluate the successes and failures of the MGGPOD package in reproducing TGRS 
data, and provide identifications for the numerous instrumental lines. 
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1. Introduction 

Due to the opacity of the Earth’s atmosphere, astronomical observations in the 7 -ray regime 
must be performed at the top of the atmosphere or in space, regions that are pervaded by intense 
radiation fields. Interactions of the particle radiations with the instrument and spacecraft materials 
result in a complex instrumental background. Discerning the much smaller signals from celestial 
sources against this strong background is the single worst problem for 7 -ray astronomy at low 
and intermediate energies of about 50 keV - 10 MeV. Thus a detailed qualitative and quantitative 
understanding of the physical processes giving rise to instrumental line and continuum backgrounds 
is crucial for most stages of 7 -ray astronomy missions, ranging from the design and development 
of new instrumentation through performance prediction and mission planning to data reduction. 

. Obtaining quantitative estimates of instrumental backgrounds is rendered v.ery difficult by the 
complexity of the physics involved. The instrumental background of a given mission has a com- 
plicated dependence on its specific radiation environment, which in turn depends on the mission’s 
orbit and epoch of operation, and on the details of the instrument design (the detector material, 
the field-of-view, active and passive shielding, and the amount and distribution of active and pas- 
sive material). For many years semi-empirical methods and scaling laws have been used to predict 
the intrumental background of a given instrument based on data from past missions and a well- 
developed understanding of the relevant physics (see e.g. Dean, Lei & Knight 1991; Gehrels 1992). 
However, the subtleties and complexities involved in scaling from one combination of radiation 
environment and instrument design to another degraded the accuracy of these predictions (Dean 
et al. 2003). 

The advent of particle transport codes capable of dependably representing the fundamen- 
tal physics involved, together with recent advances in computer processing speed, make a b initio 
Monte Carlo simulation a feasible approach for obtaining quantitative estimates of instrumental 
backgrounds (Dean et al. 2003). In this approach model spectra of the mission specific radiation 
environment are combined with a computer representation of the geometrical structure and ma- 
terial composition of the 7 -ray instrument and the rest of the spacecraft to track the trajectories 
and interactions of the incident particles and their secondaries, and to record their energy deposits 
throughout the system. One of the first Monte Carlo packages developed for this purpose was the 
University of Southampton’s GGOD suite, which provided the capability of modeling instrumental 
line and continuum backgrounds due to the delayed decay of radioactive isotopes produced in parti- 
cle interactions, and which has been applied to modeling numerous 7 -ray instruments, among them 
the GRIS balloon spectrometer, TGRS on board Wind , BATSE on board the Compton Gamma- 
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Ray Observatory , IBIS and SPI on board INTEGRAL , and BAT on board Swift (see Lei et al. 
1996; Dean et al. 2003, and references therein). Monte Carlo codes with similar capabilities have 
been developed by other groups and were validated using e.g. data from the Oriented Scintillation 
Specrometer Experiment (OSSE) on board the Compton Gamma-Ray Observatory (Dyer et al. 
1994) or from the Ge spectrometer on board the HEA 0-3 mission (Graham et al. 1997). 

Motivated by the need for accurate modeling of the instrumental line and continuum back- 
ground expected for the Ge spectrometer SPI on board the INTEGRAL mission, we have improved 
the GGOD suite and combined it with the user-friendly NASA/GSFC MGEANT package, which 
is e.g. used for SPI response simulations (Sturner et al. 2000, 2003). Both the MGEANT and 
GGOD codes are based on the widely used GEANT Detector Description and Simulation Tool 
(Version 3.21, Brian et al. 1995). In order to include photons due to the prompt de-excitation 
of excited nuclei produced in neutron captures, inelastic neutron scattering, and spallations in our 
simulations, we furthermore created the PROMPT package, which is integrated into the MGEANT 
and GGOD codes. The resulting suite of Monte Carlo packages, named MGGPOD, supports a b 
initio Monte Carlo simulations of both prompt and delayed instrumental backgrounds, including 
the plethora of instrumental 7-ray lines. The MGGPOD suite and documentation are available to 
the public for download from http : //sigma-2 . cesr .f r/spi/MGGPOD/. 

We demonstrate the capabilities of the MGGPOD suite by modeling high resolution 7-ray 
spectra recorded by the Transient Gamma-Ray Spectrometer (TGRS) on board Wind during 1995 
(earlier modeling using the GGOD suite has been reported by Lei et al. 1996; Diallo et al. 2001). 
The TGRS is a Ge spectrometer operating in the 40 keV to 8 MeV range (Owens et al. 1995). 
Due to its fine energy resolution, these spectra reveal the complex instrumental background in 
formidable detail, particularly the many prompt and delayed 7-ray lines for which we provide 
identifications. We evaluate the successes and failures of the MGGPOD package in reproducing 
the TGRS line and continuum backgrounds. A synopsis of the TRGS results presented here, and 
preliminary modeling results for the SPI spectrometer, can be found in Weidenspointner et al. 
(2004a). Recently, preliminary MGGPOD results have also been presented for the Reuven Ramaty 
High-Energy Solar Spectroscopic Imager (RHESSI, described in Smith et al. 2002) by Wunderer, 
Smith, & Weidenspointner (2004). A brief overview of the capabilites, functioning, and structure 
of the MGGPOD package has been given by Weidenspointner et al. (2004b). 

In this paper, we provide an overview on a variety of aspects pertaining to the production 
and characterization of instrumental backgrounds in 7-ray astronomy in § 2. The capabilites, 
functioning, and structure of the MGGPOD package, including the modeled physics, are described 
in § 3. Our MGGPOD modeling of TGRS data, our TGRS instrumental line identifications, and 
detailed comparisons of data and simulation are given in § 4. A summary of our results and 
concluding remarks can be found in § 5. 
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2. Instrumental Background 

In this section we provide an overview on a variety of aspects pertaining to the production and 
characterization of instrumental backgrounds in 7-ray astronomy. The instrumental background 
of a given 7-ray instrument has a complicated and complex dependence on its specific radiation 
environment and the details of the instrument design. The radiation environment, i.e. the ambient 
photon and particle radiation fields, their spectral and angular distributions, as well as their time 
histories, are a sensitive function of a mission’s orbit relative to the Earth and the geomagnetic field. 
The photons and particles that constitute the local radiation environment interact with detector 
and instrument and spacecraft structures through a wide variety of physical processes, producing 
secondary photons and particles, as well as radioactive nuclei which upon their (delayed) decay 
give rise to further secondaries. These primary and secondary photons and particles can result 
in detector triggers that pass all logical and electronical criteria required for proper events due 
to photons. from celestial sources; these triggers constitute the instrumental background. In a 
given radiation environment strength and spectral distribution of the instrumental background are 
strongly influenced by details of the instrument design such as the choice of detector material, the 
field-of-view, active and passive shielding, and the amount and distribution of active and passive 
material. 

Among the radiation fields that are potential sources of intrumental background in 7-ray 
astronomy missions are cosmic rays (Galactic cosmic rays, Solar energetic particles, and anomalous 
cosmic rays), geomagnetically trapped particles, Earth albedo radiations, diffuse cosmic X and 7 
radiation, and locally produced secondary radiation (see e.g. Stassinopoulos 1989; Klecker 1996; 
Dean et al. 2003, and references therein). An accurate description of the radiation environment 
clearly is of great importance for instrumental background estimates. 

• Galactic cosmic rays: Galactic cosmic rays primarily consist of protons (about 85% by num- 
ber), followed by a particles (^ 12%), and electrons (~ 2%); the remainder are heavier nuclei. 
For kinetic energies above about 10 GeV/nucleon solar modulation is no longer effective and 
the differential energy spectrum of cosmic rays can be approximated by an E~ 2 ‘ 7 power law. 
This particle radiation is a dominant source of instrumental background for missions that 
spend most of their time outside the Earth’s magnetosphere or at least above the radiation 
belts, but may be less significant for balloon borne experiments or missions in low-Earth or- 
bits (LEOs) because of the geomagnetic field acts as a momentum/energy filter requiring that 
a charged particle must exceed a minimum rigidity 4 * * (the so-called cut-off rigidity) to reach a 
given location within the magnetosphere. For LEOs the incident cosmic radiation therefore 
varies with time because the local cut-off rigidity changes as the spacecraft follows its trajec- 


4 The rigidity R of a particle of charge q and momentum p is given by: R—^.k cut-off rigidity of 10 GV, which 

is a “typical” value for LEOs with low inclination, implies that a cosmic-ray proton must exceed a kinetic energy of 

about 9.1 GeV to reach the spacecraft 
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tory. Outside the magnetosphere the only source of temporal variation in the incident cosmic 
radiation is solar modulation, which affects particles with kinetic energies less than a few 
GeV/nucleon. At balloon altitudes and in LEO the geomagnetic cutoff and Earth shadowing 
result in a highly anisotropic cosmic-ray intensity; outside the magnetosphere cosmic rays can 
be assumed isotropic.* 

• Solar energetic particles and anomalous cosmic rays : Solar energetic particles (SEPs) are pro- 
duced in violent energy releases on the' Sun such as flares or coronal mass ejections; during 
these times their flux can be enormously high compared to other particle radiations. SEPs 
have energies up to several 100 MeV/nucleon. Hence missons in LEOs and balloon experi- 
ments are shielded from most of them, while SEPs can be a dominant temporary source (up 
to a few days) of instrumental background for missions outside the magnetosphere. Similar to 
SEPs, anomalous cosmic rays are mainly a concern for missions outside the magnetosphere; 
because they consist mostly of heavy nuclei they can be effective producers of secondaries in 
nuclear interactions. The intensity of anomalous cosmic rays is much more strongly effected 
by the solar cycle than that of Galactic cosmic rays. 

• Geomagnetically trapped radiation: The geomagnetic field can be well approximated by a 
dipole out to a few Earth radii. Charged particles can be magnetically trapped in such a 
dipole field and stored in so-called radiation belts in the geomagnetic equatorial plane. The 
most abundant trapped particle species are protons and electrons. Protons, which are of main 
concern for instrumenal backgrounds, are trapped in a single belt, with the maximum proton 
intensity occuring at an altitude of < 1 Earth radius. The energies of the protons extend up 
to several 100 MeV; their angular distribution can be highly anisotropic (Whatts, Parnell, 
& Heckman 1989). The geomagnetic dipole is offset from the Earth’s center, therefore the 
proton belt appears to extend to lower altitudes over the coast of Brazil - this area is usually 
referred to as the South Atlantic Anomaly (SAA). Because of its intensity the trapped proton 
radiation it is avoided as much as possible by 7-ray missions. However, the SAA usually can 
not be avoided altogether for missions in LEOs. If encountered, passages through the SAA 
are a dominant source for the production of radio-nuclei in the instrument and spacecraft 
materials (Kurfess et al. 1989; Weidenspointner et al. 2001). The SAA dosage received by a 
mission in LEO strongly depends on the orbit’s altitude and inclination. 

• Earth albedo radiation : When entering the Earth’s atmosphere Galactic cosmic-ray particles 
interact violently with the air nuclei, initiating nuclear interaction cascades that ultimately 
result in the production of a multitude of secondaries of relatively low energy. The most 
important of these secondaries for 7-ray experiments at balloon altitudes and in LEO are 
photons and neutrons. At balloon altitudes, the angular distribution of these secondaries is 
highly anisotropic, and their intensity depends on the local geomagnetic cut-off rigidity and 
on the depth in the atmosphere (see e.g. Gehrels 1985, and references therein). Depending 
on energy, both albedo photons and neutrons can constitute strong and anisotropic radiation 
fields for missions in LEO, and are dominant background sources for balloon experiments. 
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The importance of the Earth’s 7 -ray albedo as a diffuse photon source increases with energy, 
above about 35 MeV it is 10-100 times more intense than the diffuse cosmic 7 radiation 
(Thompson & Simpson 1981). The energy distribution of albedo neutrons extends from 
thermal energies to several 100 MeV. The intensity of albedo 7 -rays (Thompson & Simpson 
1981; Harris, Share, & Leising 2003) and neutrons (Weidenspointner et al. 1996; Morris et al. 
1998) varies with the local geomagnetic rigidity and with the solar cycle. 

• Diffuse cosmic X and 7 radiation : The diffuse cosmic X and 7 radiation is of great astrophys- 
ical interest in its own right (e.g. Weidenspointner & Varendorff 2001), however, it constitutes 
a background against which all other observations must be made. For the purpose of back- 
ground simulations the diffuse cosmic X and 7 radiation can be considered isotropic (Kinzer 
et al. 1997; Sreekumar et al. 1998; Weidenspointner et al. 2000), and energy spectrum con- 
stant in time and unaffected by the geomagnetic field. For missions in LEO the shielding 
effect of the Earth needs to be taken into account. At balloon altitudes there is direction 
dependent attenuation by the atmosphere in addition, although the atmospheric 7 radiation 
usually dominates. 

• Internally produced secondary particles : The photons and particles of the external radiation 
fields described above interact with the instrument and spacecraft structures through a wide 
variety of physical processes (see below), producing a multitude of secondary particles such as 
photons, electrons and positrons, neutrons, protons and heavier nuclei. Because this locally 
produced secondary radiation is produced by the external radiation fields it is therefore not 
a component of the radiation environment in the same sense as the external components. 
We list here the secondary radiation as a separate component of the radiation environment 
because it is treated as such in semi-empirical calculations of instrumental backgrounds (e.g. 
Gehrels 1992). In a b initio Monte Carlo calculations, the topic of this paper, these secondary 
particles are derived from the modeling of the interactions of the external radiation fields, 
which are the only input. 

The photons and particles of the radiation environment interact through a large variety of 
processes with the instrument and spacecraft structures and contribute to the instrumental back- 
ground in different ways. Instrumental background due to external photon fields is most important 
for instruments with poor spatial resolution and a large field-of-view. The photons trigger the 
detector either by entering through the field-of-view (aperture flux) or by entering through the 
veto shield without it being triggered (shield leakage). Energetic protons (and heavier nuclei) are a 
major source of instrumental background. They may pass through the instrument and spacecraft 
suffering only slight energy losses due to ionization in the traversed materials. However, they may 
also undergo a catastrophic nuclear interaction with a single nucleus, creating a potentially large 
number of secondary particles such as pions, neutrons, protons, and light nuclei; usually a relatively 
heavy product nucleus is left behind. The secondary particles may be energetic enough to initiate 
further nuclear interactions with other nuclei, their decay (e.g. pions) can result in the generation 
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of secondary photons, electrons, and positrons, and the de-excitation of residual nuclei gives rise 
to more secondary particles (mainly photons). 

All of these processes and interactions occur on timescales that are short compared to typical 
timescales for event processing by the instrument electronics (which typically are a few to about 
100 fis) and contribute to the so-called prompt instrumental background. Depending on energy, 
veto shields can reduce prompt background components by a few orders of magnitude. Some of 
the product nuclei resulting from proton interactions may not be stable and emit further secondary 
particles when undergoing radioactive decay. Depending on the lifetime of the unstable nucleus, 
the time delay between its production and its decay may be much longer than typical instrumental 
event processing timescales, contributing to the so-called delayed instrumental background. For 
satellite missions the long-term build-up of delayed background due to long-lifed isotopes (e.g. 22 Na) 
deserves special attention (Kurfess et al. 1989; Weidenspointner et al. 2001). The best strategy to 
minimize delayed instrumental background is to minimize the amount of passive material in the 
instrument, particularly in the vicinity of the detectors. Neutrons are another major source of both 
prompt and delayed instrumental background components. Unlike protons, neutrons can traverse 
veto shielding without triggering it and produce prompt background within the detector through 
elastic and inelastic scattering (by passing energy to the recoil nucleus) as well as neutron capture 
(de-excitation of the product nucleus). Outside the veto shield the de-excitation of product nuclei 
from neutron induced nuclear interactions including neutron capture also contribute to the prompt 
background. Some of these product nuclei are unstable, and their radioactive decay adds to the 
delayed background. 

An important aspect of the instrumental background are lines. The position, strength, and 
shape of instrumental lines are relevent because they can interfere with spectroscopy of astrophysical 
7 -ray lines. Instrumental lines can result from prompt processes such as de-excitation of excited 
nuclei as well as from delayed radioactive decays. Instrumental lines can also be used as powerful 
diagnostics of the overall instrumental background, including continuum components. For example, 
7 -ray lines produced by isomeric transitions in Ge detectors can be used to estimate the neutron 
flux inside the veto shield, which then in turn allows an estimate of continuum backgrounds due to 
elastic and inelastic neutron scattering and neutron induced /3-decays (e.g. Naya et al. 1996). 
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3. The MGGPOD Monte Carlo Simulation Suite 

The MGGPOD suite is a user-friendly Monte Carlo simulation package that is applicable to 
all stages of space-based 7 -ray astronomy missions. In particular, the MGGPOD suite allows a b 
initio simulations of instrumental backgrounds - including the many 7 -ray lines - arising from 
interactions of the various radiation fields within the instrument and spacecraft materials. It is 
possible to simulate both prompt instrumental backgrounds, such as energy losses of cosmic-ray 
particles and their secondaries, as well as delayed instrumental backgrounds, which are due to the 
decay of radioactive isotopes produced in nuclear interactions. MGGPOD can also be used to study 
the response of 7 -ray instruments to astrophysical and calibration sources. The MGGPOD suite 
is therefore an ideal Monte Carlo tool for supporting most stages of 7 -ray missions, ranging from 
design, development, and performance prediction through calibration and response generation to 
data reduction; software and documentation are available to the public for download at the Centre 

f c 

d’Etude Spatiale des Rayonnements . In this publication we focus on the physics simulated by the 
MGGPOD suite. Detailed practical advice for users on how to install and use this Monte Carlo 
package, including examples, can be found in the documentation available on the MGGPOD web 
site (see also Weidenspointner et al. 2004b). 

MGGPOD is a suite of five closely integrated Monte Carlo packages, namely MGEANT, 
GCALOR, PROMPT, ORIHET, and DECAY, each of which will be described in more detail be- 
low. The MGGPOD package resulted from a combination of the NASA/GSFC MGEANT (Sturner 
et al. 2000, 2003) and the University of Southampton’s GGOD (Lei et al. 1996; Dean et al. 2003) 
Monte Carlo codes. Both were improved, and supplemented by the newly developed PROMPT 
package. The overall structure of the MGGPOD package is illustrated in Fig. 1 . Depending on 
the simulated radiation field or 7 -ray source distribution one or three steps, requiring two or three 
input files, are needed to obtain the resulting energy deposits in the detector system under study. 
In general, it is advisable to simulate each component of the radiation environment separately. MG- 
GPOD distinguishes two classes of radiation fields. Class I comprises radiation fields for which only 
prompt energy deposits are of interest, such as celestial or laboratory 7 -ray sources or cosmic-ray 
electrons. Class II comprises radiation fields for which in addition delayed energy deposits resulting 
from the activation of radioactive isotopes need to be considered. Examples for Class II fields are 
cosmic-ray protons, or geomagnetically trapped protons. 

For both of these classes, the simulation of the prompt energy deposits requires two input files: 
a mass model, and a model of the simulated radiation field. The mass model is a detailed computer 
description of the experimental set-up under study. It specifies the geometrical structure of instru- 
ment and spacecraft, the atomic and/or isotopic composition of materials, and sets parameters that 
influence the transport of particles in different materials. Each component of the radiation envi- 
ronment (and analogously for 7 -ray sources) to which the instrument is exposed is characterized 


5 http: //sigma- 2. cesr.fr/spi/MGGPOD/ 
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by three quantities: the type of the incident particles, and their spectral and angular distribu- 
tions. The prompt energy deposits are written to an output event file. For a Class I radiation field 
MGEANT is sufficient for the simulation, for a Class II radiation field an executable linking the 
MGEANT, GCALOR, and PROMPT packages is required. In case of a Class II radiation field 
there is an additional output file which lists, for all nuclei produced in hadronic interactions, the 
product nucleus 5 identity along with the geometrical mass model element in which it was produced. 
From this file isotope production rates can be computed. 

To simulate delayed energy deposits (Class II radiation field) two additonal steps need to be 
taken. These require as input the time history of the radiation field which is responsible for the 
activation, and the previously calculated isotope production rates. Based on this information first 
the activity of each isotope produced in each structural element of the mass model is determined 
by ORIHET. Then, employing the MGEANT, GCALOR, and DECAY packages, these activities 
are used to simulate the delayed energy deposits due. to radioactive decays in the instrument. 

Combining prompt and delayed energy deposits from each component of the radiation envi- 
ronment and 7-ray sources, it is possible to obtain the total energy deposited in the system as 
a function of position and time. In the following, each of the five packages that constitute the 
MGGPOD simulation suite is described. 


3.1. MGEANT 

MGEANT is a multi-purpose simulation package developed by the Low Energy Gamma Ray 
Group (LEGR) at NASA/GSFC (Sturner et al. 2000, 2003). It is based on the GEANT Detector 
Description and Simulation Tool (Version 3.21) created and supported by the Application Soft- 
ware Group, Computing Networks Division, at CERN Geneva, Switzerland (Brun et al. 1995). 
GEANT is designed to simulate the passage of elementary particles through an experimental set- 
up, which may be of considerable complexity. Although originally designed for high-energy physics 
experiments, GEANT has found applications in many other areas, including space science and 
specifically 7 -ray astronomy. Within the MGGPOD suite, MGEANT (i.e. GEANT) stores and 
transports all particles, and treats electromagnetic interactions from about 10 keV to a few TeV. 
In addition, MGEANT provides the option to use the GEANT Low-Energy Compton Scattering 
package GLECS, which provides more detailed physical models of the coherent (Rayleigh) and 
incoherent (Compton) photon scattering processes than those included in the standard GEANT 
distribution by taking into account the kinetic energy of the bound electrons (Kippen 2002). 

MGEANT was created to increase the versatility and user-friendlyness of the GEANT sim- 
ulation tool. A modular, “object oriented” approach was pursued, giving MGEANT two main 
advantages over standard GEANT. First, the instrument specific geometries and materials are pro- 
vided via input files, rather than being hard-coded. Second, several event-generation beam models 
and spectral models are available (e.g. Sturner et al. 2000, 2003). Beam models include a plane 
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wave, an isotropic radiation field, astrophysical and calibration point sources, or a user defined sky 
map; spectral models include power law and exponential spectra, line emission, or user defined 
spectra. Furthermore, MGEANT can interactively display (using CERN’s PAW++ package) the 
geometric set-up as well as the particle trajectories - a very convenient capability when creating a 
mass model or defining and verifying beam parameters. The original MGEANT offers a choice of 
output formats, in MGGPOD the only supported output file format is FITS, for which two different 
event list formats are available: in standard format the total energy deposit in each detector (and 
in the anti-coincidence system) for each event is listed; in extended format the energy deposits and 
location of each interaction, the time of interaction, and the type of the interacting particle are also 
listed. The extended format was introduced to facilitate detailed instrument response studies, such 
as the pulse-shape discrimination (PSD) system of the INTEGRAL spectrometer SPI, or event 
reconstruction algorithms for advanced Compton telescope concepts. 

MGEANT is therefore very well suited for rapid prototyping of detector systems, it can readily 
generate most of the radiation fields relevant to 7-ray astronomy, and it is set up to support detailed 
instrument response studies. The MGEANT simulation package and a user manual are available 
at the NASA/GSFC web site of the LEGR group 6 . 


3.2. GCALOR 

GCALOR (Zeitnitz & Gabriel 1994, 1999) is an interface between the CALOR89 package by 
(Gabriel et al. 1995) and the GEANT simulation tool. The CALOR89 package, designed to simulate 
calorimeter systems for high-energy physics detectors, simulates hadronic interactions down to 
1 MeV for nucleons and charged pions and down to thermal energies (10“ 5 eV) for neutrons. 
GCALOR incorporates the capability to perform hadronic interaction calculations in the GEANT 
framework by extracting the HETC (High Energy Transport Code) collision and evaporation model 
and the FLUKA model (which is already available in recent versions of GEANT) from CALOR89. 
The HETC Monte Carlo code consists of two parts: the Nucleon Meson Transport Code (NMTC) 
and the Scaling Model (which provides a smooth transition between the NMTC and FLUKA 
regimes). For low energy neutrons the MICAP (Monte Carlo Ionization Chamber Analysis Package) 
neutron code by Johnson & Gabriel (1988) has been included into GCALOR, rather than the 
M.ORSE neutron transport code utilized in CALOR89. The particle types and energy ranges 
covered by the four modules of GCALOR are: 

• NMTC: nucleons 1 MeV to 3.5 GeV, charged pions 1 MeV to 2.5 GeV 

• Scaling Model: nucleons and charged pions 3 GeV to 10 GeV 


6 http://lheawww.gsfc.nasa . gov/do cs/gam cosray/legr/mgeant/mgean t.html 
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• FLUKA: nucleons and charged pions above 10 GeV, and for all energies for particle types 
not implemented in CALOR 

• MICAP: neutrons 10 -5 eV to 20 MeV 

Hadronic interactions between particles that are heavier than individual nucleons, such as deuterium 
or helium nuclei, cannot be simulated with GCALOR. However, continuous energy losses, for 
example due to ionization of the matter traversed, are taken into account for all charged particles. 
The GCALOR package is publicly available at the Universitat Mainz 7 . 

When combining GCALOR with GEANT (or MGEANT), all particles are stored and trans- 
ported by GEANT. Electromagnetic interactions are simulated by GEANT, hadronic interactions 
are simulated by GCALOR. The GEANT plus GCALOR package therefore extends the capabilities 
of the standard GEANT tool to include hadronic interactions of charged particles down to 1 MeV 
and of neutrons down to thermal energies. Equally important, this package provides access to the 
energy deposits from all interactions as well as to isotope production anywhere in the geometrical 
set-up. 

Originally, GCALOR/MICAP utilized only ENDF/B (Evaluated Nuclear Data Files version 
B) neutron data. Unfortunately, these data, available e.g. at the National Nuclear Data Center 
(NNDC) at Brookhaven National Laboratory 8 , do not cover all elements or isotopes. In particular, 
individual neutron cross-sections for the five natural Ge isotopes were missing. These cross-sections, 
which are clearly of great importance for simulations of instrumental backgrounds in Ge spectrome- 
ters, were generated based on JENDL (Japanese Evaluated Nuclear Data Library 9 ) data by Zeitnitz 
(2001, priv. comm.). 

Some of the low-mass residual nuclei produced in spallations as modelled by GCALOR are Very 
neutron rich, particularly if the incident projectile is energetic (e.g. a primary cosmic-ray proton). 
Some of these neutron rich nuclei may be spurious, some have been identified in nuclear experiments 
and found to have short lifetimes of less than 1 s. The neutron rich nuclei are unstable against the 
emission of one or several neutrons, a process that affects the number and energy distribution of 
secondary neutrons. To include this source of secondary neutrons into our simulations, MGGPOD 
checks whether light product nuclei ( Z < 8) are unstable against neutron emission, exploiting 
existing data whenever available (Firestone 1996), and mandates neutron emission if appropriate, 
converting the unstable original product nucleus into the highest-mass isotope that has no neutron 
exit channel. 


7 http: / /wswww. physik.um-mainz.de / zeitmtz/gcalor/gcalor.html 

8 http://www.nndc.bnl.gov/ 

9 http://wwwndc. tokai.jaerLgo.jp/jendl/jendlhtml 
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3.3. PROMPT 

The PROMPT package was created for simulating the prompt de-excitation of excited nuclei 
produced by neutron capture, inelastic neutron scattering, and spallation 10 . The PROMPT pack- 
age consists of a data base containing information on the de-excitation of a large number of nuclei, 
and of code to access these data and to generate random samples of de-excitation particles. When 
simulating prompt hadron-induced backgrounds with MGGPOD, the MGEANT, GCALOR, and 
PROMPT packages are linked into a single executable (see Fig. 1). Each time a secondary nu- 
cleus is produced in a nuclear interaction as simulated by GCALOR, PROMPT is called to model 
the de-excitation cascade. In GCALOR, prompt photons are generated for a few hadronic inter- 
actions, however, in general these photons are “statistical” and continuous - they do not reflect 
the actual discrete nuclear levels and the well-know selection rules for transitions between them. 
This approximate treatment proved sufficient for modeling calorimeters for high-energy physics 
detectors (Zeitnitz & Gabriel 1994), but is clearly insufficient for detailed modeling of instrumental 
7 -ray backgrounds. Therefore any prompt photons generated by GCALOR are replaced with those 
returned by the PROMPT package. 

We developed the PROMPT package in the spirit of the compound nucleus model, in which 
it is assumed that the reaction process can be separated into independent incoming and outgoing 
channels with a well-defined intermediate compound state. We rely on the GCALOR package for a 
complete specification of the incoming channel in the three classes of prompt reaction, i.e. neutron 
capture, inelastic neutron scattering, and spallation. Beyond a determination of the daughter 
nucleus the outgoing channel is, however, at best incompletely specified by GCALOR; in particular, 
de-excitation photons are either missing or incorrect, as described above. In other words, GCALOR 
generates a known daughter nucleus in an unknown excited state with unknown decay properties. 
In addition, the distribution of excited states left behind by the three types of prompt reaction is 
not generally known - unlike the case of radioactive decay, where the excited state in the daughter 
has a well-known probability from extensive measurements of branching ratios (see § 3.5). To create 
the PROMPT package we therefore had to derive (simple) recipes for selecting an excited state 
in the daughter nucleus and then modeling its de-excitation into known lower energy levels. Our 
approach is to use a probabilistic recipe to specify the intermediate excited state of the compound 
daughter nucleus, and a further ansatz to specify the de-excitation of the excited state into known 
lower levels. When generating de-excitation particles (usually photons), we do not take into account 
Doppler broadening due to the motion of the nucleus during de-excitation, which is relevant for 
short-lived levels with T\/ 2 < 0.5 ps (Evans et al. 2002). 


10 Although spallation is physically defined in terms of the evaporation model (see § 3.3.2), for the purpose of 
MGGPOD it is defined to comprise all hadronic interactions other than neutron capture which result in a product 
nucleus that is different from the target nucleus. We will, however, assume the validity of the evaporation model for 
all such interactions in our treatment 
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3.3.1. Neutron Capture and Inelastic Neutron Scattering 

In the case of neutron capture and inelastic neutron scattering, the excitation energies involved 
are fairly low and the processes involved are fairly well understood. We selected the excited state 
of the compound nucleus at random from a level density formula, by which the probability of 
obtaining an energy level depends on its spin J, parity tt and excitation energy U. The value 
of U was obtained from GCALOR from the kinematics of the incoming channel reaction. The 
distribution of spins and parities follows the angular part of the Bethe level density formula as 
parametrized by Mughabghab and Dunford (1998); this angular part contains a weak dependence 
on U which we neglected (the overall density contains an exponential in y/U , but we are only 
interested in the relative probabilities of different J and 7T values). Having selected an excited state 
in this way, we assumed that it coupled to known states of lower excitation energy (from the ENSDF 
database 11 ) by an electric dipole (El) transition - the spins and parities of the lower states were 
required to obey the well-known El selection rules for coupling to J 7 ’’, and if none were permitted 
we assumed magnetic dipole (Ml) transitions. The probability of a transition from an excited state 
to a permitted level x was assumed to be proportional to ({/ - U x ) 3 ) which holds approximately 
for dipole transitions (Blatt & Weisskopf 1952) - we assume that the evaporation model is valid 
for spallation reactions as defined above. The specific prescriptions for generating prompt photons 
for these two neutron-induced reactions are outlined in the following, and illustrated in schematic 
form in Figs. 2 and 3. 

For neutron captures on nuclei with A > 3 (M)GEANT and GCALOR specify the kinetic 
energy of the incident neutron, the target and the product nucleus, and the recoil energy E r of 
the product nucleus. We estimated the excitation energy U of the product nucleus as the sum 
of its neutron separation energy S n (taken from Firestone 1996, typically several MeV), which 
is the excitation energy after capturing a thermal neutron, and the kinetic energy of the incident 
neutron Ek ; the recoil energy of the product nucleus typically is only a few eV and therefore 
negligible. We obtain the probability of a compound nucleus level having given values of from 
the Mughabghab & Dunford (1998) level density formula as described above; however, we make 
the approximation that it is always evaluated at excitation energy U ~ 5 n , which is justified by the 
weak variation of the J * distribution on U . The de-excitation cascade which follows was treated 
by two different ansatze y depending on the value of J T \ those values which are compatible with 
thermal neutron capture were simply assigned the branching ratios from the extensive thermal 
neutron database maintained by NNDC, while all other values were assumed to de-excite by 
the El mechanism described above. Thermal neutrons couple to the target nucleus by low partial 
waves in the incoming channel (we assumed s-wave), so that the corresponding J n in the compound 
nucleus are easily specified, given the target nucleus’ ground state; for example, in the common 
case of an even-even target (0 + ), an s-wave corresponds to in the compound nucleus. Our 
motivation was to exploit the very large amount of thermal neutron experimental data available 


11 Evaluated Nuclear Structure Data Files (ENSDF) are e.g. available at the NNDC 
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in this special case. The product nuclei for which the prompt de-excitation cascade after neutron 
capture is modeled in these ways are listed in Table 1. 

For inelastic neutron scattering on nuclei with A > 4 (M)GEANT and GCALOR specify, 
among other quantities, the kinetic energy of the incident neutron E' K , the recoil nucleus and its 
kinetic energy E ri as well as that of the scattered neutron E f k . The excitation energy of the recoil 
nucleus U is assumed to be U ~ Ek — E k — E r . In general, the GCALOR cross-sections are accurate 
enough to reflect the lowest levels. However, above about 1-2 MeV the derived excitation energies 
approach a “continuum” . Modeling of the de-excitation cascade depends on the excitation energy 
U - If it is less than or equal the energy of the highest known level with transition data in ENSDF, 
Eh, then the initial level of the de-excitation cascade is chosen to be the level with transition data 
whose energy is closest to U. If U is greater than Eh then we assume a single initial transition from 
U down to any of those known levels for which spin-parity and transition probability data exist. 
The probability for this initial transition is calculated from the standard spin-parity distribution 
and the dipole selection rules as indicated above. The spin-parity distribution is evaluated at a 
reference energy E re j chosen for each nucleus to be approximately where the inelastic neutron 
scattering cross-section peaks according to JENDL, which characteristically is around 10 MeV. As 
in the case of neutron capture (where we used U ~ 5 n ), we assume that the calculated probability 
distribution for the initial transition from E rt j applies to all excitation energies U greater than Eh- 
This is justified by the weak dependence of this distribution on U. In the de-excitation cascade 
from this level at excitation energy U the energy of the initial transition is increased by U - E re f- 
The recoil nuclei for which the prompt de-excitation cascade after inelastic neutron scattering is 
modeled are listed in Table 1. 

The existence of asymmetrically broadened triangular or “sawtooth” shaped features extending 
in energy above the lines from the lowest energy neutron inelastic scatters has been recognized by 
several authors in both laboratory (Bunting k Kraushaar 1974) and flight (Wheaton et al. 1989; 
Evans et al. 2002; Weidenspointner et al. 2003) data, and are seen in our data also. They arise from 
the summed energy depositions by the recoiling Ge nucleus and by the de-excitation photon (s). 
These features are automatically predicted by (M)GEANT and GCALOR and furnish a comparison 
between the simulation and the behavior of the spectrometer. 


3.3.2. Spallation 

For spallation reactions (M)GEANT and GCALOR specify the type and kinetic energy and 
momentum of the incident particle, the target nucleus, and the types, energies and momenta of all 
reaction products. However, because GCALOR conserves energy and momentum for any interaction 
only on average, but not rigorously for an individual case, it is not possible to obtain a meaningful 
estimate of the excitation energy U of any reaction product from kinematics. Nor are experimental 
cross-sections available for production of such individual levels, which are extremely numerous. 
However, some genera] rules about the distribution of U values follow from the evaporation model, 
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which is the generalization of the compound nucleus model to particle emission (Weisskopf & 
Ewing 1940). Proton-induced spallation reactions at high energy leave the compound nucleus in a 
very broad distribution of excitation energies, which can be thought of as a thermal distribution 
~ exp(~U/kT) for some temperature parameter T (unlike heavy-ion collisions, which tend to 
select collective states such as the high-spin yrast levels: Galin 2001). If T exceeds the nucleon 
separation energy one or more such nucleons are likely to be emitted (neutrons above S n being 
favored, since there is no Coulomb barrier). The generation of a daughter nucleus by GCALOR 
implies that neutron emission has ceased, which in turn implies that T < S n . Characteristic values 
of S n are several MeV; the fact that the thermal distribution has a substantial "tail” at energies 
above T implies that neutron emission only ceases when T falls far enough below S n for the ” tail” 
to become negligible. Experiments attempting to achieve high T in residual nuclei left behind by 
proton reactions on nuclei between Ag and Au have attained T ~ 5 MeV, with values 3-4 MeV being 
characteristic (Ledoux et al. 1998). On this basis we assumed that the known excited states U x in 
spallation product nuclei were distributed according to the formula (2 J x + l)exp(-f7 x /3 MeV). A 
sketch of the procedure followed once particle emission has stopped is given in Fig. 4. 

A key assumption behind this type of argument is that nuclear energy levels can be treated 
statistically, which is recognised to become less and less valid as the nuclear mass A falls below 
about 30 (Holmes et al. 1976). We found that, in practice, this thermal approximation cannot 
be used at all for A < 20 since far too many lines are predicted in the simulation. On the other 
hand, the relatively small number of individual levels involved means that data on experimental 
cross-sections to individual levels of A < 20 nuclei are in some cases quite complete. Ramaty, 
Kozlovsky, & Lingenfelter (1979) reviewed these cases, and we have obtained from their compilation 
the distribution of excited states for the nuclei in Table 2 when these are produced by spallation. 
The de-excitation of all other nuclei with A < 20 was omitted from the simulation. 

The statistical model is also expected to break down for low incident particle energies (below 
a few MeV). In particular, this applies to spallations that are induced by secondary neutrons and 
protons. In these cases individual levels, rather than a statistical distribution parameterized by T, 
will be excited. In addition, the difference in level distributions produced by protons and neutrons 
(which we have otherwise ignored) will be important. 


3.4. ORIHET 

The original ORIHET program was developed for the GGOD Monte Carlo suite (Lei et al. 
1996; Dean et al. 2003). It is an adaptation of the Oak Ridge Isotope Generation and Depletion code 
ORIGEN, which was designed to calculate the build-up and decay of activity in any system for which 
the nuclide production rates are known. To do so, ORIHET has a built-in data base containing 
the half-lives and decay channels of a large number of radio-isotopes. Currently, ORIHET returns 
the activity at a given time for two different radiation histories: constant irradiation, and “cooling 
down” after constant irradiation. In the first case, it is assumed that nuclei are produced at a 
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constant rate for a given time period. This model is useful for simulating instrumental backgrounds 
for missions in radiation environments that are relatively stable, or for deriving average instrumental 
backgrounds for missions exposed to variable radiation fields as are encountered e.g. in LEOs. In 
the second case, it is assumed that nuclei are produced at a constant rate for a given period of 
time during which production and decay compete, then the production ends and the radio-isotopes 
decay for a given “cool-down” period. This model is e.g. useful for estimating the instrumental 
background due to transits through the SAA. 

The original ORIHET radionuclide data base did not include all isotopes relevant for simula- 
tions of instrumental backgrounds in 7 -ray experiments. Specifically, for some isotopes which can 
be produced in their ground state as well as an isomeric state, such as 24 Na, the isomeric state 
was absent from the original data base. For MGGPOD, we have changed ORIHET to include all 
isomers missing in the original data base which produce significant 7 - ray lines identified in spectra 
of TGRS or other Ge spectrometers. In addition, we corrected the treatment of a-decays. 


3.5. DECAY 

As ORIHET, the original DECAY package was created for the GGOD suite (Lei et al. 1996; 
Dean et al. 2003). The DECAY package is the analog of PROMPT for simulating radioactive 
decays. It consists of a data base containing information on the decays of a large number of 
radioactive isotopes, and of code to access these data and to generate random samples of decay 
particles. Linking the MGEANT, GCALOR, and DECAY packages into a single executable (see 
Fig. 1) enables one to simulate the decays of a given radio-isotope with a given activity over a given 
time period at any location in the geometric set-up. 

To generate the DECAY data base, information concerning the decay schemes, such as branch- 
ing ratios for and electron capture (EC) decays, internal transition (IT), and a-decay chan- 

nels, as well as the energy levels and 7 -ray to internal conversion electron branching ratios, was 
taken from ENSDF. These data have been supplemented with information on X-ray and Auger 
electron fluorescence yields and energies as given in Firestone (1996). 

In developing MGGPOD, we improved and expanded the DECAY package. The number of 
isotopes in the data base was increased, and data base and code upgraded to simulate isomeric 
levels properly (these sometimes give rise to significant instrumental lines in 7 -ray detectors). 
Information on the emission of internal conversion electrons is very incomplete in ENSDF data. 
For a few isotopes relevant for Ge detectors we corrected the data files (e.g. for the 691 keV E0 
transition in 72 Ge, for which photon emission is not possible). 
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4. Modeling Data of the Ge Spectrometer TGRS 
4.1. Instrument and Radiation Environment 

The Transient Gamma- Ray Spectrometer (TGRS) on board the Wind spacecraft was primarily 
designed to perform high- resolution spectroscopy of transient 7-ray events, such as 7-ray bursts or 
solar flares (Owens et al. 1995). The detector itself consists of a 215 cm 3 high purity n-type Ge 
crystal sensitive to energies in the 20 keV - 8 MeV band, kept at its operating temperature of 
85 K by a passive radiative cooler constructed mainly of Be and Mg. Some shielding in the soft 
X-ray range, mainly against intense solar flare X-rays, is provided by a 30 mil (0.762 mm) thick 
sheet of Be-Cu alloy around the sides of the cooler. The TGRS detector is located on the south- 
facing surface of the rotating cylindrical Wind body, which points permanently toward the southern 
ecliptic pole. The spectrometer has no active shielding and is permanently exposed to ^ 1.87T sr of 
the southern hemisphere which is unobstructed by the cooler.' A 1 cm thick Pb occulter attached 
to the Wind body exploits the spacecraft rotation in order to modulate the signal from the ecliptic 
plane; as seen from the detector it occults a band 90° long by 16° wide sweeping out the whole 
ecliptic in Wind’s 3 s rotation period. The chemical elements mentioned, together with Al which 
was used for most structural components, are expected to be among the main sources of prompt 
and delayed 7-ray line backgrounds. 

Since its launch on November 1, 1994, Wind has been following unusual and highly elliptical 
orbits; halo orbits around the Earth-Sun Lagrangian LI point in the first half of the mission, 
trans-lunar Earth orbits in the second half (Acuna et al. 1995). Hence Wind spent virtually 
the whole mission in interplanetary space, well away from near-Earth radiation backgrounds such 
as geomagnetically trapped particles and Earth albedo radiations. The radiation environment 
experienced by TGRS therefore is dominated by two components: diffuse cosmic hard X and 7 
radiation, and Galactic cosmic rays, unmodulated by the Earth’s magnetic field. Both radiation 
fields can be considered isotropic. The diffuse cosmic hard X and 7 radiation is contant in time, 
significant secular changes of Galactic cosmic rays due to solar modulation occur on time scales of 
a few months. These changes were observed by TGRS in several 7-ray lines (Harris et al. 2001). 
There have been very few interruptions in the data stream, generally caused by brief passages of 
the trapped radiation belts around perigee or by triggering of a special mode of data collection by 
solar flares and 7-ray bursts. 

The background spectra recorded by TGRS provide an ideal test for the MGGPOD package. 
The (relative) simplicity of the instrument design facilitated the development of an accurate model 
for Monte Carlo simulations (this also applies to the Wind spacecraft). The stability and isotropy 
of the radiation environment, which can be represented by onfy two dominant components, also 
greatly simplify Monte Carlo simulations, and at the same time render quantitative comparisons of 
TGRS spectra with simulations more rigorous. Finally, the fine energy resolution of the detector 
reveals the plethora of instrumental lines in great detail, providing us with sensitive tests on the 
numerous interaction channels through which Galactic cosmic rays can deposit their energy in the 
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instrument and spacecraft structures. The resolution of TGRS at 500 keV was nominally about 
3 keV full width at half maximum (FWHM), which was achieved in the early months of the mission 
(Harris et al. 1998); thereafter resolution degraded due to accumulated damage from cosmic-ray 
impacts (Kurczynski et al. 1999). The line profiles also became distorted, with marked tails on the 
low-energy wings, and energy calibration became more difficult due to gain shifts. We limit our 
analysis to the period 1995 Jan.-May, when these problems were negligible, except for small low- 
energy tails particularly on the highest energy lines (see § 4.3). The TGRS background spectra were 
binned in 1 keV energy channels over the entire 20 keV - 8 MeV range. However, due to problems 
with saturation of electronics components after very high energy deposits by heavy Galactic cosmic 
rays the energy range of TGRS for scientific studies had to be limited to 40 keV - 8 MeV. In 
addition, the 210-260 keV range is contaminated by electronic artefacts, which appear as broad 
features in the spectrum. Our reference Jan.-May, 1995 TGRS spectrum is shown as part of Fig. 5. 

The detector cannot distinguish individual energy deposits which occur within the peaking 
time of the electronic, which is about 5 /xs. In general, nuclear transitions within the Ge crystal 
from excited states to the ground state will be detected at the sum of the energies if they proceed 
through one or more intermediate levels. An exception occurs for transitions via states with lifetimes 
r > 5 //s, which may survive to be detected separately; we refer to these states as isomers for the 
purpose of working with TGRS. This is particularly relevant for the de-excitation of 73m Ge: the de- 
excitation of the 66.7 keV level (r — 0.72 s) proceeds via a level at 13.3 keV with r = 4.3 fis. If the 
13.3 keV level takes long to de-excite, the decay will be registered as two distinct events, giving rise 
to lines at 13.3 keV and 53.4 keV; if it de-excites quickly the two photons will be summed giving rise 
to a single line at 66.7 keV; for intermediate cases the energies are only partially summed and the 
event will fill in the interpeak region between 53.4 keV and 66.7 keV. Radioactive j3~ or /3 + decays 
within the detector do not give rise to line features in the spectrum; the energy of the /3-particle 
is distributed continuously and summed with the coincident 7 rays of discrete energy. Electron 
capture /?-decays within the Ge crystal do produce multiple line features. Electron captures are 
accompanied by X-rays and/or Auger electrons arising from the filling of a vacancy in one of the 
atomic sub-shells; lines will appear at the sum of the nuclear transition energy and the binding 
energy of the atomic sub-shells. Certain lines will appear in the TGRS spectrum which are not 
due to cosmic-ray effects, notably a strong line at 1460 keV from a 40 K source which was flown for 
calibration purposes, and lines from the decay chains of U and Th occurring naturally in traces, 
e.g. in the Be cooler and shield. 


4.2. Models and Simulation 

In modeling the TGRS background by Monte Carlo simulation approximations regarding the 
mass model representing the instrument and the Wind spacecraft, as well as the components of 
the radiation environment and their time history, are inevitable. All elements of the TGRS instru- 
ment described above and of the Wind spacecraft are sources of prompt and delayed instrumental 
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backgrounds. The required mass model has to give a faithful representation of the geometrical 
structure and its atomic/isotopic compositon, particularly for the Ge detector and its vicinity. 
With increasing distance from the detector the details of the mass distribution are less critical; 
however, we were careful to conserve the total mass of the spacecraft components and to represent 
their atomic/isotopic composition, as these are important factors for the generation of secondary 
particles and residual nuclei in hadronic interactions. Our model of the TGRS detector is based 
on the mass model created by Seifert et al. (1995) for response matrix generation. This early mass 
model replicated well the geometrical structure and the masses of the TGRS instrument above and 
around the Ge crystal as described by Owens et al. (1995). However, it did hot specify the atomic 
and isotopic compositions of materials in sufficient detail for our purpose, it did not include instru- 
ment components below the Ge crystal, such as electronics boxes, and it only provided the crudest 
representation of the spacecraft. We supplemented and improved the early TGRS mass model by 
referring to original technical drawings and other documentation, and by close examination of spare 
parts. Assuming that the latest documents on the detector mass that were available to us are the 
most accurate, we can account for about 96% of the mass of the total instrument (18.1 kg out of 
18.9 kg). However, the “mass deficit” in and close to the detector, i.e. the areas which are most 
important for the simulation, we can account for only 86% of the material (6 kg out of 7 kg). Our 
approach to creating a model for the Wind spacecraft and other scientific instruments was similar. 
For example, the basic geometry and structure of the spacecraft body are described in Harten & 
Clark (1995). We can account for about 86% of the documented spacecraft mass in our mass model. 

The radiation environment of TGRS was approximated by its two dominant components: 
diffuse cosmic hard X and 7 radiation, and cosmic radiation. Both were assumed to be isotropic 
and constant in time. We used an analytic approximation by Gruber et al. (1999) to model the 
spectrum of the diffuse cosmic photon radiation in the 30 keV to 9 MeV energy range. Due to 
limitations of GCALOR (see § 3.2) the only hadronic particle component of the cosmic radiation 
that could be simulated were protons (the contribution of cosmic-ray electrons to TGRS spectra 
was found to be negligible). We used a proton spectrum (energy range 10 MeV to 226 GeV) that 
was calculated for early 1995 by Moskalenko (2002, priv. comm.) with the GALPROP Galactic 
cosmic-ray progagation code, taking into account solar modulation using a steady-state drift model 
(Moskalenko et al. 2002). This spectrum was assumed to be representative of the Jan.-May 1995 
time period, during which the TGRS data we are modeling were recorded. During this time period 
the activity of some isotopes produced by cosmic radiation built up, and hence their contribution 
to the TGRS data increased with time. However, the TGRS spectrum we are modeling represents 
an average over this time period. In our simulations of the instrumental background due to delayed 
radioactive decays we approximate this time average with a “snap shot” in March 1995 after four 
months of cosmic-ray irradiation, i.e. we use isotope activities that are calculated assuming four 
months of constant irradiation with our model cosmic-ray proton spectrum. 
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4.3. Identifications of TGRS Instrumental Lines 

One of the main goals for developing the MGGPOD suite was to create a Monte Carlo tool 
capable of modeling the many instrumental lines present in 7-ray detectors. To assess the extent to 
which this goal was reached for TGRS we attempted to identify its more than 200 observable lines 
and spectral features. The TGRS spectra were analyzed using the GASPAN 12 gamma spectrum 
analysis program. GASPAN allowed us to characterize the shape of the instrumental lines as a 
function of energy, following parameterizations given by Phillips & Marlow' (1976). The photopeak 
component of each line was described by a Gaussian, the low-energy “wing”, whose importance 
increased with energy, was modeled by the sum of two exponential tails representing effects from 
multiple Compton scatterings and incomplete charge collection. The line width and tail parameters 
were determined as a function of energy using the strongest and cleanest lines. This energy de- 
pendent line shape was then employed to analyse TGRS spectra, particularly to resolve the many 
complex regions resulting from. closely spaced or blended lines. When analyzing intrinsically broad 
lines arising from the prompt de-excitation of very short-lived levels (T x / 2 < 0-5 ps, Evans et al. 
2002) allowance had to be made for their exceptional width, which is caused by Doppler broadening 
due to the motion of the nucleus during de-excitation. 

The results of our spectral analysis are summarized in Tables 3 and 4, in which the lines are 
identified and the count rates in the observed spectrum are presented. Lines or blended features for 
which the GASPAN fit determined a significance < 5a are not included. Most of the strong lines had 
already been identified with nuclear transitions in previous work with Ge spectrometers (Wheaton 
et al. 1989; Bartlett 1994; Evans et al. 2002; Weidenspointner et al. 2003), with which our results 
are in general agreement. However, the level of detail with which our simulation was performed 
allowed us to suggest that in many cases multiple transitions can be identified contributing to 
previously undetected blends. Thus in Tables 3 and 4 two or more transitions are frequently 
assigned to an observed line feature. The order in which the transitions are listed corresponds 
to their ranking as contributors to the line strength, according to the simulation, with transitions 
whose contribution was < 10% according to the simulation being omitted. Table 5 lists identified 
triangular or “sawtooth” shaped features from inelastic neutron scattering off Ge atoms in the 
detector crystal, as described in § 3.3.1. 

Narrow (unresolved) lines were assigned to Table 3 if one or more identification could be 
made, from earlier work or from the simulation. If no certain identification could be made they 
were put in Table 4. Lines which were visibly broader than the TGRS energy resolution were 
treated as blends, keeping in mind that an intrinsically broad component might also contribute. 
If multiple identifications were possible, sufficient (in theory) to explain the line width, the line 
was included in Table 3. Two indications were used to assign a blend to the "uncertain identity” 
category (Table 4). First, a blend which was obviously broad enough to contain multiple lines might 


12 The software, developed by F. Riess, and documentation are available under http://ftp.leo.org/download/ 
pub/science/physics/software/gaspan/; alternatively, F. Riess can be contacted directly (friedrich@die-riessens.de) 
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correspond to only one line in the simulation (or none). Second, a comparison of the observed and 
predicted line strengths might show the simulation seriously failing to reproduce an observed line. 
Although there is considerable scatter in the reliability of the simulations even for well-identified 
narrow lines (see below), we somewhat arbitrarily assigned blends where the ratio of simulated to 
observed line strength fell below 25% to Table 4. 

There are some obvious deficiencies in this qualitative procedure by which Tables 3 and 4 
were compiled. There are considerable quantitative uncertainties in the simulations, as we shall 
see, preventing the ranking of transitions and the 10% cutoff from being totally reliable. The 
simulation may even omit lines altogether. Such cases belong in Table 4, but may appear in 
Table 3 if the line appears narrow or if > 25% of its observed strength is predicted. Several cases 
of this kind may be indicated by measurements of the line energy in Table 3 which are inconsistent 
with the transition energies. However, these cases may also be due to problems inherent in the use 
of the rather complicated function (Phillips & Marlow 1976) for fitting the line profiles. 


4.4. Comparison of Simulation and Data 

A comparison of the Jan -May 1995 TGR.S spectrum with a MGGPOD simulation is shown in 
Fig. 5. The simulated instrumental background components are: prompt background due to cosmic- 
vzy proton interactions and prompt de-excitations (green), delayed backgrounds from radioactive 
decays in the TGRS instrument and the Wind spacecraft (purple), and background due to diffuse 
cosmic X and 7 rays (blue). The sum of all simulations is depicted in red, the data are shown in 
black (the broad features in the 210-260 keV region are electronic artefacts). 

4-4-1 . Continuum and Inelastic Neutron Scattering Features 

The MGGPOD simulation reproduces very well the overall shape and magnitude of the actual 
background. This is illustrated in Fig. 6, which depicts the ratio between the simulation and the 
measured spectrum. Disregarding line features, this ratio exhibits very little trend with energy and 
has a value of about 0.85. The only exception occurs at the lowest energies, below about 200 keV, 
where the ratio increases with decreasing energy up to a value of 1.2. 

Below about 200 keV diffuse cosmic photons are found to be the dominant source of instru- 
mental background in TGRS, as expected for a wide field-of-view instrument. Prompt cosmic-ray 
proton induced events are the main background component at energies above about 200 keV, and 
practically the sole background above about 4 MeV. The relative dominance of the prompt back- 
ground in TGRS reflects the lack of any veto shielding. The prompt background comprises different 
processes, among them direct ionization in the detector by cosmic-ray protons and energy losses by 
secondaries, including energy deposits due to elastic and inelastic scattering of secondary neutrons 
off Ge nuclei in the detector crystal. The prompt background is not a pure continuum, but features 



- 22 - 


numerous lines; these will be addressed in § 4.4.2. 

As discussed in § 3.3.1, inelastic neutron scattering gives rise to characteristic “saw-tooth” or 
triangular-shaped broad features in the spectrum. Their absolute and relative strengths, and less 
so their detailed shape, are sensitive to the flux and spectral distribution of secondary neutrons in 
the Ge crystal at energies around 1 MeV. In Table 5 we attempt to make comparisons between our 
measurements of the triangular inelastic neutron scattering features and the simulations. Fig. 5 
makes it clear that upon these features there are superposed a large number of strong lines, that 
underlying them there must be a continuum of somewhat uncertain shape, and that the features 
sometimes overlie each other. Nevertheless, there is sufficient resemblance between the simulated 
and observed recoil features for comparisons to be made - lower and higher energy bounds and 
lines can be consistently identified, hence simulation and observation can be treated consistently. 

Our method involved simply identifying the lower and upper energy bounds of the recoil fea- 
tures, summing the counts between these limits, and subtracting the estimated counts contributed 
by the lines and the continuum. The uncertainties in the values in Table 5 are overwhelmingly 
dominated by the uncertainties in the true values of line and continuum count rate. We quote 
ranges of count rate in each recoil feature corresponding to the most extreme estimates of the line 
and continuum strengths. The extreme low continuum estimate assumes a fiat continuum having 
the flux corresponding to the upper energy bound; the high continuum estimate is a linear inter- 
polation between the lower and higher energy bounds. The line measurements within the recoil 
feature energy ranges (Table 3) were assumed to be upper limits because of the possibility of con- 
tamination by the recoil feature. Of course, upper line and continuum limits correspond to lower 
limits on the strength of the recoil feature in Table 5, and vice versa. 

In Table 5, note that the upper and lower energy bounds are somewhat arbitrarily chosen from 
one feature to the next, due to the interference of strong lines and other recoil features. Thus the 
results should not be compared from one feature to the next - the important point is that theory 
and observation can be compared within the same feature, since they are treated consistently. 

We are left with the conclusion, based on the figures in Table 5 (where comparable) that 
the simulation and the measurement in general agree within 30% - an exception is the one above 
834 keV, which suffers exceptionally strong interference from the lines at 844 and 847 keV. This 
estimate is supported by the lack of any discontinuities in the ratio of simulation and actual data 
at the recoil features’ energy ranges in Fig. 6 . There is a tendency for the theory to over-predict 
the strength of the recoil features. This might result from systematic errors in any of the various 
quantities that are relevant for the simulation of the recoil features, such as the flux of secondary 
neutrons - which depends e.g. on the accuracy of the mass model and of hadronic cross-sections in 
general, and the Ge neutron cross-sections in particular. 

Delayed radioactive decays in the TGRS instrument and the Wind spacecraft also give rise to 
continuum background ( 7 -ray lines from radioactive decays are discussed below). The dominant 
contributor are /3-decays in the Ge crystal. As described in § 4.1, /3“ or decays within the 
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Ge detector result in a pure continuum background (except if isomeric levels are involved in the 
de-excitation of the daughter nucleus) because the continuously distributed energy of the /3-particle 
is summed with the coincident 7 rays of discrete energy. In the case of TGRS and Wind , continuum 
background due to radioactive decays in the spacecraft is relatively less important because m/r 2 , 
with m being a mass element and r the distance from the detector, is lower for the spacecraft 
structure than for instrument parts. The delayed background exceeds the diffuse cosmic photon 
background above about 400 keV, and cuts off at about 4 MeV. 


4 -4 Lines 

The simulation is also very successful in modeling the more than 200 lines that are observed 
in the TGRS spectrum. Most (about 87%) of the lines are reproduced, with the ratio of simulated 
and actual line count rates clustering around a value of 1 with no trend in energy. The simulation 
produces a few spurious lines, and fails to reproduce a few lines that are present in the data. 

A quantitative comparison of simulated and observed line strengths is shown in Fig. 7, in the 
form of a scatter plot of the ratio between them as a function of line energy. The distribution of 
the numerical values of these ratios are shown in Fig. 8. Overall, the observed lines are reproduced 
within a factor ^ 2.5. There is no trend with energy in Fig. 7, and the mean ratio is very close 
to 1, suggesting that our methods of line simulation (§§ 3.5, 3.3) are sound to a first level of 
approximation. There are enough lines (205) for the results to be broken down into subsets in 
search of interesting systematic deviations from these overall conclusions. 

If the simulated-to-observed line strength ratios are broken down by reaction type it is found 
that lines from radioactive decays are on average well reproduced by the simulation, while spallation 
reaction line strengths are somewhat overestimated (Fig. 8, dotted line). Leaving aside possible 
problems of isotope production from the incoming channel cross-sections in GCALOR, this effect 
results solely from a systematic overestimate of the strengths of lines from low-mass (20 < A < 30) 
nuclei, by a factor of nearly 2 (Fig. 8, dashed line). The selected spallation product nuclei with 
A < 20 which were included in the simulation (Table 2) were not overproduced 
!!! TBD: 12C a problem ???? update scatter plot !!!. 

This indicates that the thermal approximation for spallation product de-excitations (§ 3.3) becomes 
increasingly inadequate with decreasing mass number below A = 30, as might be expected from 
nuclear level statistics. Spallation reactions overall do not show a wider spread than the others (a 
factor ^ 2.5 in Fig. 8), nor do they show any trend with line energy. This indicates that, to the 
extent that the thermal approximation is valid, the choice of temperature 3 MeV is justified. 
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5. Summary and Conclusion 

We have described the capabilites, functioning, and structure of the MGGPOD Monte Carlo 
simulation suite, which is based on the widely used GEANT Detector Description and Simulation 
Tool (Version 3.21, Brun et al. 1995). MGGPOD is a user-friendly simulation package that is 
applicable to all stages of space-based 7-ray astronomy missions. In particular, the MGGPOD 
package allows a b initio simulations of instrumental backgrounds, both continuum and 7-ray lines, 
that arise from interactions of the various radiation fields within the instrument and spacecraft 
materials. It is possible to simulate both prompt instrumental backgrounds, such as energy losses 
of cosmic-ray particles and their secondaries, as well as delayed instrumental backgrounds, which are 
due to the decay of radioactive isotopes produced in nuclear interactions. The MGGPOD package 
and documentation are available to the public for download at the Centre d’Etude Spatiale des 
Rayonnements. 

To demonstrate its capabilities, we employed the MGGPOD suite for modeling high resolution 
7-ray spectra recorded by the Transient Gamma-Ray Spectrometer (TGRS) on board Wind during 
1995. We found that both the continuum and the 7-ray line background are very well reproduced. 
Regarding the continuum background, the ratio between simulation and flight data is about 0.85, 
except for the lowest energies below about 200 keV, where the ratio increases with decreasing 
energy up to a value of 1.2. The simulation even reproduces well asymmetric, “sawtooth” or 
triangular-shaped broad spectral features that are due to inelastic neutron scattering off Ge nuclei 
in the detector. Regarding background lines, the simulation reproduces about 87% of the observed 
background lines in TGRS, with the ratio of simulated and actual line count rates clustering around 
a value of 1 with no trend in energy. Our results for continuum and lines are as good as or better 
than those obtained by Dyer et al. (1994) for OSSE, by Dyer et al. (1998) for Mars Observer, or 
Graham et al. (1997) for HEAO-3 using other Monte Carlo tools. 

Below about 200 keV, diffuse cosmic photons are found to be the dominant source of instrumen- 
tal background in TGRS. The simulation overestimates this component, with the largest deviation 
of 20% occurring at the lowest energies. We attribute the overestimate mainly to uncertainties in 
the mass distribution around the Ge detectors. Although we made a substantial effort to describe 
the detector as accurately as possible in our mass model, we may have missed material in and close 
to the detector (as described in § 4.2). Small changes in the amount and distribution of mass close 
to the detector have a large effect on the simulated count rate at a few tens of keV because of 
the strong photoelectric absorption at these energies, and also affect the probability of Compton 
scattering of higher-energy photons into the detector. It is also possible that the employed model 
spectrum of the diffuse cosmic photons (Gruber et al. 1999) is too high. . . 

??? few words on uncertainty in XRB - find recent XRB review ??? 

The short-fall of the simulation above 200 keV is expected. Currently, GCALOR is not ca- 
pable of simulating hadronic interactions of a particles, hence the contributions of cosmic-ray a 
particles to the prompt and the delayed background could not be modeled. A rough estimate of 
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this contribution can be obtained with an empirical formula for the ratio of the total a to pro- 
ton cross-sections by Ferrando et al. (1988).- This ratio does not vary much with the mass of the 
target nucleus between Al and Ge, and has an average value of about 1.7. By number, the ratio 
of cosmic-ray as to protons is about 0.14; hence we can estimate that if cosmic-ray a particles 
could be included in the simulation, they would increase the prompt and delayed background by 
roughly 24% - enough to account for the current short-fall by about 15%. The effect of a induced 
spallations on the activation of specific isotope species is very difficult to estimate, since the ratio 
of a to proton cross-sections for specific spallation products strongly varies with energy (Ferrando 
et al. 1988). Similarly, when modeling data from the OSSE instrument on board the Compton 
Gamma-Ray Observatory , Dyer et al. (1994) concluded that a particles account for about 20% of 
the calculated activation due to cosmic rays. In addition, there are deficiencies in the mass model, 
since the total mass accounted for in the instrument and the spacecraft falls short of the docu- 
mented mass budget by about 14% (§ 4.2). Less mass will result in fewer interactions and hence 
less prompt background as well as less activation and therefore less delayed background. Except 
for the lowest energies, shielding of the detector from 7 -rays originating in different parts of the 
detector and spacecraft is not expected to be significantly affected by the moderate mass deficit. 

The simulation reproduces about 87% of the observed background lines in TGRS, with the 
ratio of simulated and actual line count rates clustering around a value of 1 with no trend in 
energy. On average, lines from radioactive decays are well reproduced. We see a small trend of 
overproduction for lines from de-excitation of spallation products with decreasing mass. This trend 
is expected given the simple statistical method that is used for generating de-excitation photons. 

The simulation of lines from de-excitations involving isomeric levels is influenced by the choice 
for the detector time constant (§ 4.1). In our simulations, we set this parameter to Td = 5 ps\ hence 
the results for lines from the isomeric levels of 73 Ge (the 13.3 keV level has a lifetime of 4.3 ps) are 
particularly sensitive to this parameter. As can be seen in Figs. 5 and 6 , both lines at 53.4 keV 
and 66.7 keV are overestimated, with the sum peak being relatively too strong. In the simulation, 
partial summation of energy deposits by the electronics, which fills the interpeak region, is not 
modeled, thereby overestimating the sum peak. However, to reproduce the shape of the 73m Ge 
complex better in the simulation, the parameter Td needs to be set to a value smaller than 5 ps. 

There are a few lines which allow us to assess the overall quality of the simulation. The 511 keV 
line, which results from a multitude of processes (e.g. decay of 1 r + produced in hadronic interactions, 
electromagnetic showers, radioactive decays) is reproduced within a factor of 0.74, indicating that 
there are no severe problems with modeling the diverse physics involved. The 2.223 MeV line 
from thermal neutron capture on H, the features from inelastic neutron scattering off Ge, and the 
Doppler broadened line at 2.211 MeV line from inelastic neutron scattering on 27 A1 are reproduced 
within factors of 1.4, < 1.3, and 1 . 1 , respectively. This indicates that flux and energy spectrum of 
secondary neutrons throughout TGRS and Wind are fairly well reproduced, which in turn implies 
that there are no severe problems with the treatment of neutron production, propagation, and 
interaction. 
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The remarkable performance of MGGPOD in modeling TGRS data encouraged us and others 
to employ this Monte Carlo suite to other instruments. Preliminary results for the spectrometer SPI 
on board INTEGRAL have been reported by Weidenspointner et al. (2004a). MGGPOD reproduces 
well early-mission SPI data; however, remaining discrepancies imply that in the case of this much 
heavier and complex mission (compared to TGRS /Wind) there are deficiencies in the production 
and/or thermalization and capture of secondary neutrons. This has lead us to an ongoing systematic 
evaluation of the neutron cross-sections in GCALOR, whose results will be presented in a future 
publication of a re-visited SPI background simulation. Recently, MGGPOD has also been applied 
to RHESSI by Wunderer, Smith, & Weidenspointner (2004). These results are still preliminary, 
but provide us with clues for improving the modeling of the time-dependent activation during 
SAA passages. Two of us (G. W. and M. J. H.) will use MGGPOD to estimate instrumental 
background and sensitivity of the space based gamma- ray lens MAX (von Ballmoos et al. 2004). 
The performance and user-friendliness of MGGPOD also attracted the various instrument teams 
studying concepts for an advanced Compton telescope; they selected MGGPOD as their baseline 
simulation tool for predicting the performance and sensitivity of the various designs 
!!! TBD reference if possible !!!. 


We are grateful to David Palmer, Helmut Seifert, and Juan Naya for helpful discussions about 
the TGRS instrument. We are indepted to Prof. Friedrich Riess for enhancing GASPAN according 
to our needs, and for user support. 
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Fig. 1. — A flow chart illustrating the overall structure of the MGGPOD Monte Carlo simulation 
suite. The various simulation packages (shown in boxes) and input and output files (shown in 
ellipses and round-edged boxes) are explained in the text 

Fig. 2. — Schematic representation of our treatment of (n, 7 ) reactions. The incoming neutron (left, 
bold arrow) is captured in a level at excitation energy S n + E * in the compound nucleus, where 
S n is the neutron separation energy and Ek the neutron’s kinetic energy. Spin a and parity n are 
assigned to this level using the level density formula of Mughabghab & Dunford (1998) to obtain a 
random value (inset (a)). In the general case the level is then allowed to decay electromagnetically 
to known levels of the nucleus (arrows down) which have compatible spin and parity according to 
the electric dipole selection rules; the probability of transition to any level is proportional to E*. In 
the special case (b) where the J * of the excited level are compatible with s-wave neutron capture 
the downward transition probabilities are obtained directly from measured thermal neutron capture 
data (ENSDF).. 

Fig. 3. — Schematic representation of our treatment of neutron inelastic scattering. The excitation 
produced in the target nucleus is known from the GCALOR particle tracking. Two cases arise. 

(a) The incoming neutron (left, bold arrow) excites a level in the target which lies above the 
highest known level E^. The excited level is assigned a spin and parity at random from the level 
density distribution (curve inset, top right). The level density formula is evaluated at energy 
E re f ~ 10 MeV, approximately where the inelastic cross section peaks. A de-excitation transition 
is assigned to the highest known level with J * compatible with an El transition (dashed arrow). 

(b) The incoming neutron deposits energy of a known level (or between known levels); the nearest 
level in energy is assumed to be excited. In either case, once a known level has been occupied the 
cascade of transitions down to the ground state is taken from the ENSDF database. 

Fig. 4. — Schematic representation of our treatment of spallation reactions. The incoming particle 
(left, bold arrow) at very high energies leaves the compound nucleus in a thermal distribution 
of excited states which loses energy by particle emission. When it has lost enough energy that 
only a small ’’tail” remains above the neutron separation energy S n , particle emission ceases. The 
thermal distribution function shown at right (rotated) then applies (the ’’tail” being shaded), and 
the known levels (full) in the remnant nucleus are populated according to a temperature T = 3 MeV. 
De-excitations from these known levels (an example at energy U x is shown) are then taken from 
the ENSDF nuclear database. We do not consider the unknown levels (dash) because there are no 
data for them in ENSDF. 

Fig. 5. — A comparison of the Jan .-May 1995 TGRS spectrum with a MGGPOD simulation. 
Details are given in the text. The broad features in the data in 210-260 keV are electronic artefacts 

Fig. 6. — The ratio of the MGGPOD simulation and the Jan.-May 1995 TGRS spectrum 

Fig. 7. — Ratio between simulated and observed line strengths for TGRS background spectrum 
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accumulated during 1995 Jan.-May. Top arrows - lines predicted by simulation but not observed. 
Bottom arrows - observed lines without simulated equivalents (mostly from A < 20 spallation 
product nuclei not included in simulation). The lines are well reproduced, typically to within a 
factor of 2.5 

Fig. 8. — Distribution of ratios of simulated to observed line strengths. Full histogram - all lines. 
Dotted histogram - de-excitation lines from spallation reactions. Dashed histogram - de-excitation 
lines from spallation reactions yielding nuclei with mass 20 < A < 30. The symbol “>>” indicates 
the bin containing lines predicted but not observed (infinite ratio) 
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Table 1. Isotopes for which prompt de-excitation photons are generated after neutron capture 

and/or inelastic neutron scattering. 


Process Product Isotopes 


Neutron Capture (71,7) 10 Be, 28 Al, 49 Ti, 53 Cr, 55,57 Fe, 64,66 Cu, 71 *73,74,75,77 Q e 

Inelastic Neutron Scattering (n, n f 7) 27 A1, 54 ’ 56 ^ 7 Fe 3 63,65 Cu, 70 ^2i73,74,76Q e 


Table 2. Nuclei produced by spallation for which de-excitation line probabilities were obtained 
from Ramaty, Kozlovsky, & Lingenfelter (1979) 


Process Product Isotopes 

Spallation (A < 20) 10 B, n ’ 12 C, 14 N, 15 ’ 16 0, 19 F 
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Table 3. Identified lines; blends of identified lines 


TGRS line Count rate Transition ID Lab Comment 4 

energy* 3 s _1 (levels) b,c energy* 3 


53.4 a,s ~ 15% buildup on 73 As t x / 2 . 
53.9 a From 65 Ga(/3 + ). 

53.0 s 


53.6 

0.71 

73m Ge(67-13) 
65m Zn(54-g.s.) 
58 Co(53-g.s.) 
46 Sc(52-g.s.) 

59.2 

0.100 

60,11 Co(59-g.s.) 

? 74m Ga(60-g.s.) 

66.5 

0.65 

73ra Ge(67-13-g.s.) 

88.0 

0.25 

69m Ge(87-g.s.) 

24 Na(563-472) 

92.6 

0.88 

67m Zn(93-g.s.) 
67m Zn(93-g.s.) + L 

119.5 

0.051 

72m Ga(120-16-g.s.) 

139.7 

0.47 

75m Ge(140-g.s.) 

? 57 Fe(136-14-g.s.) + L 

143.4 

0.13 

57 Fe(136-14-g.s.) + K 
46ra Sc(143-g.s.) 

159.3 

0.055 

47 Ti(159-g.s.) 

68 As(158-g.s.) 

175.7 

0.031 

71 Ge(175-g.s.) + L 
71 Ge(175-g.s.) 

185.7 

0.26 

67 Zn(185-g.s.) + L 
71 Ge(175-g.s.) + K 

194.2 

0.24 

67 Zn(185-g.s.) + K 

198.2 

1.47 

71m Ge(198-175-g.s.) 

19 F(197-g.s.) 

239.1 

0.165 

212 Bi(239-g.s.) 

271.2 

0.053 

44m Sc(271-g.s.) 
228 Th (328-58) 

278.0 

0.026 

208 Pb(3475-3198) 

228 Th(1153-874) 

284.2 

0.0023 

61 Ni(283-g.s.) + L 


52.0 

s 


58.6 

s 


59.7 

s 


66.7 

a } $ 

~ 5% buildup on 73 As t x i 2 . 
Two-step transition. 

86.8 

a, s 

From 69 As (/?+) 

91.0 

s 


93.3 

a,s 

From 67 Ga(EC). 

94.5 

a 

From 67 Ga(EC). 

119.5 

(2,5 

Two-step transition. 

139.7 

a, 5 

From 75 Ga(/3“). 

137.3 

a 

From 57 Co(EC). 

143.6 

a 

60% buildup on 57 Co ry 2 . 

142.5 

s 


159.4 

a, 5 

From 47 Sc(/?~). 

158.1 

5 


176.4 

a, 5 

From 7l As(EC). 

174.9 

n 


185.8 

a 

From 67 Ga(EC). 

186.1 

a 

From 71 As(EC). 

194.2 

a 

From 67 Ga(EC). 

198.4 

a, 5 

From 71 As(/3 + ). Two-step transition. 

197.1 

5, a 

~ 2% buildup follows solar modulation 

238.6 

r 

232 Th series. 

271.1 

s 


270.2 

r 

232 Th series. 

277.4 

r 

232 Th series. 

279.0 

r 

232 Th series. 

284.0 

a 

From 61 Cu(EC). 


Table 3 — Continued 


TGRS line' 
energy b 

Count rate 
s” 1 

Transition ID 
(levels) b,c 

Lab 

energy b 


Comment 3, 

291.2 

0.0046 

61 Ni(283-g.s.) + K 

291.3 

a 

From 61 Cu(EC). 

296.1 

0.0153 

214 Bi(295-g.s) 

295.2 

r 

238 U series 

301.8 

0.0082 

67 Zn (394-93) + L 

301.4 

a 

From 67 Ga(EC). 

304.5 

0.015 

75m As(304-g.s.) 

303.9 

s 


309.6 

0.062 

67 Zn (394-93) + K 

309.9 

a 

From 67 Ga(EC). 

320.4 

0.0056 

51 V(320-g.s.) 

320.1 

a,s 

From 51 Ti(y0"), 51 Cr(EC). 



228 Th(1154-832) 

321.6 

r 

232 Th series 

325.6 

0:0103 

51 V(320-g.s.) + K 

325.6 

a 

From 51 Cr(EC). 

329.2 

0.0056 

228 Th(328-g.s.) 

328.0 

r 

232 Th series. 



■ ? 69 Ga(3l9-g.s.) + K 

329.0 

a 

From 69 Ge(EC). 

331.9 

0.0076 

21 Na(332-gs) 

331.9 

s 


338.5 

0.013 

59 Ni(339-g.s) 

339.4 

s 




228 Th(396-58) 

338.3 

r 

232 Th series. 

350.7 

0.088 

21 Ne(350-g.s.) 

350.7 

s,a 

From 21 Na(/? + ). 



214 Bi(352-g.s.) 

351.9 

r 

238 U series. 

373.2 

0.0066 

43 Ca(373-g.s.) 

372.8 

a 

From 43 K (/?“), 43 Sc(EC). 

377.1 

0.0024 

52m Mn(378-g.s.) 

377.7 

s 


389.8 

0.0110 

2S Mg(975-585) 

389.7 

s,a 

From 25 Na {fir). 



214 Po(1764-1378) 

387.0 

r 

238 U series. 



214 Po(21 19-1730) 

389.1 

r 

238 U series. 

394.2 

0.0077 

67 Zn(394-g.s.) + L 

394.7 

a 

From 67 Ga(EC), 

402.9 

0.044 

67 Zn(394-g.s.) + K 

403.2 

a 

From 67 Ga(EC). 

416.8 

0.026 

26 Al(417-g.s.) 

416.9 

s 


426.9 

0.0034 

73m As(428-g.s.) 

428.3 

s 


439.2 

0.21 

23 Na(440-g.s.) 

440.0 


From 23 Ne(^-), 23 Mg(/3+). 



69m Zn(439-g.s.) 

438.6 

s 


450.8 

0.021 

23 Mg(451-g.s.) 

450.7 

s 




25 Al(452-g.s.) 

451.5 

s 


472.2 

0.056 

24m Na(472-g.s.) 

472.2 

s 


477.1 

0.022 

7 Li(478-g.s.) 

477.6 

s 


511.0 

2.26 

e + annihilation 

511.0 


~ 2% buildup follows solar modulation. 

538.0 

0.0042 

59 Fe(1750-1211) 

537.4 

s 
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Table 3 — Continued 


TGRS line 
energy 1 * 

Count rate 
s -1 

Transition ID 
(levels) b,c 

Lab 

energy b 


Comment a 

564.8 

0.0089 

228 Th(1531-969) 

562.5 

r 

232 Th series 



S4 Cr (3786-3222) 

563.7 

a 

From 54 V (/?"). 

575.1 

0.021 

69 Ga(574-g.s.) + L 

575.4 

a 

From 69 Ge(EC). 

584.0 

0.125 

69 Ga(574-g.s.) + K 

584.5 

a 

From 69 Ge(EC). 



25 Mg(585-g.s.) 

585.0 

s 




22 Na(583-g.s.) 

583.0 

s,a From 22 Mg(/3 + ). 



208 Pb(3198-2615) 

583.2 

r 

232 Th series. 

608.0 

0.030 

74 Ge(596-g.s.) + K 

606.9 

a 

From 74 As(EC). 



214 Po(609-g.s.) 

609.3 

r 

238 U series. 

656.7 

0.0067 

20 F(656-g.s.) 

656.0 

s 




61 Ni(656-g.s.) + L 

657.0 

a 

From 6l Cu(EC). 



61 Ni(656-g.s.) 

656.0 

s 


664.2 

0.0070 

61 Ni(656-g.s.) + K 

664.3 

a 

From 61 Cu(EC). 



214 Po(1275-609) 

665.5 

r 

238 U series 

670.4 

0.0103 

63 Cu(670-g.s.) 

669.6 

s 




38m Cl(671-g.s.) 

671.4 

s 


718.4 

0.025 

10 B(718-g.s.) 

718.3 

s 


721.4 

0.0111 

46 Sc(774-52) 

721.9 

s 


727.2 

0.0101 

212 Po(727-g.s.) 

727.3 

r 

232 Th series. 



228 Th(1123-396) 

726.9 

r 

232 Th series. 

744.0 

0.0069 

52 Cr(3114-2370) 

744.4 

a 

From 52 Mn(EC). 



234 U (786-43) 

742.8 

r 

238 U series. 

749.9 

0.0060 

56 Co(1720-970) 

750.0 

a 

From 56 Ni(EC). 

767.0 

0.0115 

234 U(810-43) 

766.4 

r 

238 U series. 



214 Po(1378-609) 

768.4 

r 

238 U series. 

783.3 

0.0041 

24 Na(1345-563) 

781.4 

s 




50 Cr(783-g.s.) 

783.3 

s 


794.5 

0.0036 

228 Th(1123-328) 

794.9 

r 

232 Th series. 



27 A1(3004-2211) 

793.0 

s 


803.4 

0.0153 

2O6 Pb(803-g.s.) 

803.1 

s 


811.2 

0.0176 

58 Fe(811-g.s.) + L 

811.6 

a 

50% buildup on 58 Co 



58 Fe(811-g.s.) 

810.8 

a,s 

From 58 Co(EC). 
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Table 3 — Continued 


TGRS line 
energy^ 

Count rate 
s -1 

Transition ED 
(levels) b,<: 

Lab 

energy* 5 


Comment a 

818.2 

0.055 

58 Fe(811-g.s.) + K 

817.9 

a 

66% buildup on 58 Co r^. 

829.9 

0.011 

26 Al(1058-228) 

829.4 

s 


835.3 

0.014 

54 Cr(835-g.s.) + L 

835.5 

a 

From 54 Mn(EC). 



54 Cr(835-g.s.) 

834.8 

a 

From 54 V(/T). 



228 Th(1023-187) 

835.7 

r 

232 Th series. 

842.7 

0.080 

27 Al(844-g.s.) 

843.7 

s, a 

From 27 Mg(/?“). 



54 Cr(835-g.s.) + K 

840.8 

a 

From 54 Mn(EC). 

847.0 

0.036 

56 Fe(847-g.s.) 

846.8 

s, a 

23% buildup on S6 Co r^. 

854.0 

0.0093 

56 Fe(847-g.s.) + K 

853.9 

a ' 

~ 25% buildup on 56 Co 

859.3 

0.0068 

208 Pb(3475-2615) 

860.6 

a 

232 Th series. 

871.1 

0.017 

24 Na(1341-472) 

869.2 

s 




24 Na(1346-472) 

874.4 

$. 




69 Ga(872-g.s.) + L 

873.3 

a 

From 69 Ge(EC). 

882.7 

0.038 

69 Ga(872-g.s.) + K 

882.3 

a 

From 69 Ge(EC). 



? 21 Na(3680-2798) 

881 

s 


890.1 

0.027 

22 Na(891-g.s.) 

890.9 

s 

' 



46 Ti(889-g.s.) 

889.3 

a 

10% buildup on 46 Sc 

899.0 

0.0092 

204 Pb(899-g.s.) 

899.2 

s 


911.3 

0.0138 

228 Th (969-58) 

911.2 

r 

232 Th series. 



61 Ni(909-g.s.) 

908.6 

a 

From 61 Co(/?“) 

931.1 

0.0041 

55 Fe(931-g.s.) 

931.3 

s 


936.5 

0.0186 

214 Po(1543-609) 

934.1 

r 

232 Th series. 



52 Cr(2370-1434) 

935.5 

a 

From 52 Mn(EC). 



? 18 F(937-g.s.) 

937.2 

s 


942.3 

0.0046 

52 Cr(2370-1434) + K 

941.5 

a 

From 52 Mn(EC). 

955.1 

0.0020 

27 Mg(1940-985) 

955.3 

s 


962.0 

0.017 

63 Cu(962-g.s.) 

962.1 

s 




228 Th(1023-58) 

964.8 

r 

232 Th series. 

968.8 

0.0075 

228 Th(969-g.s.) 

969.0 

r 

232 Th series. 

974.5 

0.0121 

25 Mg(975-g.s.) 

974.4 

s, a 

From 25 Na(/?~). 

983.4 

0.043 

27 Mg(985-g.s.) 

984.6 

s 




48 Ti(984-g.s.) 

983.5 

3, a 

From 48 V(/3 + ). 
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Table 3 — Continued 


TGRS line 
energy* 3 

Count rate 
s~ x 

Transition ID 
(levels) b ’ c 

Lab 

energy* 3 


Comment* 1 

988.9 

0.0157 

48 Ti(984-g.s.) + K 

988.5 

a 

From 48 V (EC). 



25 Mg(1965-975) 

989.9 

s 


1001.2 

0.0182 

234 U(1045-43) 

1001.0 

r 

23S U series. 



? 7O Ga(1003-g.s.) 

1002.6 

s 


1006.0 

0.0033 

26 A1 (3074-2070) 

1004.1 

s 




53 Cr(1006-g.s.) 

1006.1 

s 


1014.1 

0.059 

27 Al(1014-g.s.) 

1014.4 

5,a 

From 27 Mg (/?") 



26 A1(2069-1058) 

1011.7 

s 


1021.6 

0.021 

10 B(1740-718) 

1021.7 

$, a 

From 10 C (/?+). 

1025.3 

0.0030 

2oo Pb(1027-g.s.) 

1026.5 

s ) a 

From 200 Bi(EC). 

1049.3 

0.020 

66 Zn(1039-g.s.) + K 

1048.9 

a 

From 66 Ga(EC). 

1056.8 

0.00164 

2o F(1057-g.s.) 

1056.8 

a 

From 20 O(/3“). 

1063.4 

0.0079 

2o7m Pb(1064-g.s.) 

1063.7 

s 


1077.9 

0.0037 

68 Zn(1077-g.s.) + L 

1078.6 

a 

From 68 Ga(EC). 

1086.8 

0.0056 

68 Zn(1077-g.s.) + K 

1087.1 

a 

From 68 Ga(EC). 

1107.9 

0.0128 

69 Ga(1107-g.s.) + L 

1108.1 

a 

From 69 Ge(EC). 



7l Ge(1096-g.s.) + K 

1106.6 

a 

From 71 As(EC). 

1117.2 

0.086 

69 Ga(1107-g.s.) + K 

1117.1 

a 

From 69 Ge(EC). 



65 Cu(1116-g.s.) + L 

1116.6 

a 

From 65 Zn(EC). 



65 Cu(1116-g.s.) 

1115.5 

s 


1124.6 

0.041 

65 Cu(1116-g.s.) + K 

1124.5 

a 

45% buildup on 65 Zn Tj/ 2 . 



21 Ne(2867-1746) 

1121 

s 


1130.4 

0.020 

26 Mg(2938-1809) 

1129.7 

s 


1157.0 

0.0128 

44 Ca(1157-g.s.) 

1157.0 

a 

From 44 Sc(y3 + ). 



214 Po(2699-1543) 

1155.6 

r 

238 U series. 

1173.2 

0.0126 

60 Ni(1173-g.s.) 

1172.9 

s 

From 60 Co 

1227.2 

0.0021 

42 Ca(2752-1525) 

1227.7 

a 

From 42 Sc(/? + ). 

1237.5 

0.011 

214 Po(1847-609) 

1238.1 

r 

238 U series. 



56 Fe(2085-847) 

1238.3 

s 


1245.8 

0.00078 

56 Fe (2085-847) + K 

1245.5 

a 

From 56 Co(EC). 



14 C(7341-6094) 

1248 

s 


1274.2 

0.027 

22 Ne(1275-g.s.) 

1274.5 

s y a 

From 22 Na(/3 + ). 
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Table 3 — Continued 


TGRS line 
energy b 

Count rate 
s- 1 

Transition ID 
(levels) b ’ c 

Lab 

energy b 


Comment 3, 

1284.1 

0.0023 

47 Ti( 1444-159) 

1284.9 

s 


1303.7 

0.00130 

47 Ti (2749-1444) 

1304.6 

s 


1311.7 

0.013 

48 Ti(2296-984) 

1312.1 

a 

From 48 V (/?+). 

1316.4 

0.0064 

48 Ti(2296-984) + K 

1317.1 

a 

From 48 V(EC). 

1326.6 

0.0049 

63 Cu(1327-g.s.) 

1327.0 

s 


1332.3 

0:0112 

60 Ni(1333-g.s.) 

1332.5 

s 

From 60 Co(P~). 

1346.5 

0.0087 

69 Ga(1337-g.s.) + K 

1347.0 

a 

From 69 Ge(EC). 



24 Na(1345-g.s.) 

1344.7 

s 




64 Ni(1346-g.s.) 

1345.8 

s 


1368.5 

0.126 

24 Mg(1369-g.s.) 

1368.6 

a, 5 

From 24 Na(/T). 

1393.4 

0.00166 

21 Ne(1746-351) 

1396 

s 


1399.7 

0.0042 

22 Na(1984-583) 

1401 

s 


1407.5 

0.0070 

214 Po(2017-609) 

1408.0 

r 

238 U series. 

1412.8 

0.0049 

63 Cu(1412-g.s.) 

1412.1 

s 


1433.1 

0.0141 

52 Cr(1434-g.s.) 

1434.1 


From 52 Mn(EC). 



234 U(1435-g.s.) 

1435.4 

r 

238 U series. 



? 27 A1 (4410-2982) 

1428.1 

s 


1460.4 

0.040 

40 Ar(1461-g.s.) 

1460.8 


40 K calibration source. 

1480.4 

0.0028 

65 Cu(1482-g.s.) 

1481.8 

s 


1508.9 

0.0024 

24 Na(1512-g.s.) 

1512.3 

s 




214 Po(2119-609) 

1509.2 

r 

238 U series. 

1525.7 

0.0086 

22 Na(1528-g.s.) 

1528.1 

s 




? 42 Ca(1525-g.s.) 

1524.7 

a 

From 42 Sc(/3 + ). 

1546.9 

0.00174 

63 Cu(1547-g.s.) 

1547.0 

s 


1553.0 

0.0025 

50 Ti(1554-g.s.) 

1553.8 

a 

From 50 V(EC). 

1571.4 

0.0016 

2O1 Pb(2507-990) 

1570.8 

s 


1649.6 

0.0018 

26 A1 (2069-417) 

1652 

s 

2 lines 

1697.0 

0.0022 

27 Mg(1698-g.s.) 

1697.9 

s 


1721.5 

0.0142 

27 A1 (2735-10 14) 

1720.3 

s 

Broad 

1726.4 

0.0075 

214 Po(1730-g.s.) 

1729.6 

r 

23S U series. 



24 Na double esc. 

1732.0 

s 


1758.4 

0.0062 

214 Po(1764-g.s.) 

1764.5 

r 

238 U series. 
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Table 3 — Continued 


TGRS line 
energy b 

Count rate 
s'* 1 

Transition ID 
(levels) b,c 

Lab 

energy b 


Comment 3, 



25 Mg(2738-975) 

1762.9 

s 


1778.5 

0.0061 

28 Si(1779-g.s.) 

1778.9 

a 

From 28 Al(/3~). 

1792.4 

0.00104 

2S Mg(3405-1612) 

1793.4 

s 


1809.0 • 

0.049 

26 Mg(1809-g.s.) 

1808.6 

5, a 

From 26 Na(/?~). 



56 Fe(2658-847) 

1810.8 

a 

From 56 Mn(/3 _ ). 

1844.8 

0.00176 

214 Po(1847-g.s.) 

1847.4 

r 

238 U series. 

1861.0 

0.00132 

63 Cu(186l-g.s.) 

1861.3 

s 


1871.2 

0.00031 

24 Na(3217-1345) 

1872.0 

s 


1882.9 

0.0016 

1S N(7155-5271) 

1884.8 

s 


1902.2 

0.0017 

24 Na(3745-1846) 

1899.0 

s 




69 Ga(1892-g.s.) + K 

1901.8 

a 

From 69 Ge(EC). 

1947.6 

0.00115 

69 Ga(1924-g.s.) + K 

1934.4 

a 

From 69 Ge(EC). 

1966.0 

0.00115 

25 Mg(1965-g.s.) 

1964.5 

s 


1981.6 

0.0056 

18 0(1982-g.s.) 

1982.0 

s 


2032.3 

0.0015 

69 Ga(2024-g.s.) + K 

2034.0 

a 

From 69 Ge(EC). 

2055.1 

0.0012 

25 Mg(5462-3405) 

2056.4 

s 




23 Mg(2051-g.s.) 

2051 

s 


2081.3 

0.0014 

22 Ne(3357-1275) 

2082.5 

5, a 

From 22 F(/?“). 

2129.5 

0.0022 

11 B(2124.7-g.s.) 

2124.5 

s 




26 Mg(3941-1809) 

2132.0 

s 


2166.8 

0.0029 

22 Ne(5523-3357) 

2165.9 

s,a 

From 22 F(/3~). 



38 Ar(2167-g.s.) 

2167.4 

a 

From 38 K (/?+). 

2183.5 

0.00085 

6 Li(2186-g.s.) 

2186 

s 




17 O(3055-871) 

2184.5 

a 

From 17 N(/3~). 

2209.5 

0.0275 

27 Al(2211-g.s.) 

2211.0 

s 

Broad 



214 Po(2204-g.s.) 

2204.2 

r 

238 U series. 

2220.3 

0.0171 

p(n,i)d direct capture 

2223 

n 


2241.3 

0.0051 

24 Mg single esc. 

2243.0 

a 

From 24 Na (/?“). 

2312.0 

0.0030 

14 N(2313-g.s.) 

2312.6 

s 


2436.4 

0.00077 

21 Ne(2789-351) 

2438 

s 


2508.0 

0.0109 

26 Mg(4318-1809) 

2509.6 

s 


2614.9 

0.0139 

208 Pb(2614.6-g.s.) 

2614.5 

r 

232 Th series. 
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Table 3 — Continued 


TGRS line 
energy b 

Count rate 
s ' 1 

Transition ID 
(levels) b ’ c 

Lab 

energy b 


Comment a 



20 Ne(4248-1634) 

2613.8 

s 


2753.7 

0.028 

24 Mg(4123-1369) 

2754.0 

a,s 

From 24 Na(/?“). 

2761.4 

0.0045 

66 Zn(2752-g.s.) + K 

2761.6 

a 

From 66 Ga(EC). 

2981.8 

0.0038 

27 Al(2982-g.s.) 

2981.8 

s 

Broad. 



23 Na(2982-g.s.) 

2981.9 

s 


3003.1 

0.0120 

27 Al(3004-g.s.) 

3004.0 

s 

Broad. 

3336.0 

0.00175 

13 C single esc. 

3342.2 

s 


3404 

0.009 

12 C double esc. 

3416.0 

s 

Broad. 

3683.9 

0.0015 

13 C(3684-g.s.) 

3683.9 

s 


3800.4 

0.0029 

66 Zn(3792-g.s.) + K 

3801.2 

a 

From 66 Ga(EC). 

3853.5 

0.00171 

13 C(3854-g.s.) 

3853.2 

s 




? 24 Mg(5235-1369) 

3866.2 

s 


3920 

0.012 

12 C single esc. 

3927.0 

s 

Broad. 

4090.0 

0.0017 

66 Zn(4086)-g.s.) + K 

4095.9 

a 

From 66 Ga(EC). 

4122.3 

0.00094 

24 Mg(4122.6-g.s.) 

4122.7 

s,a 

Sum peak. From 24 Na(/3“). 

4246.8 

0.0017 

15 N double esc. 

4247.2 

s 


4306.3 

0.0013 

66 Zn(4296-g.s.) + K 

4305.5 

a 

Both from 66 Ga(EC). 



66 Zn(4296-g.s.) + L 

4297.1 

a 


4431 

0.015 

12 C(4439-g.s.) 

4438.0 

s 

Broad. 

4588.9 

0.0015 

14 N single esc. 

4593.8 

s 


4757.3 

0.0017 

1S N single esc. 

4758.2 

s 


5106.0 

0.0037 

16 0 double esc. 

5106.6 

s 




14 N(5106-g.s.) 

5104.9 

s 


526913 

0.0026 

15 N(5270-g.s.) 

5269.2 

s 


5289. 

0.0021 

15 N(5298-g.s.) 

5297.8 

s 

Broad ??? 

5298.0 

0.00048 

ls N(5298-g.s.) 

5297.8 

a 

From 15 C(/T) 

5420. 

0.0013 

? 

? 

? 

Broad 

5617.8 

0.0044 

16 0 single esc. 

5617.6 

s 


6128.7 

0.0052 

16 0(6129-g.s.) 

6128.6 

s 


6242.1 

0.00006 

? 

? 

? 


6723.6 

0.00005 

? 

? 

? 



a Key to reaction types: a — /?-decay after activation; s — spallation followed by prompt 
de-excitation; n — (n, 7 ) or {n,n!) followed by prompt de-excitation; r — natural radioactivity. 

b All energies in keV 

c K and L are atomic sub-shell binding energies in cases of electron capture (EC). 
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Table 4. Blends containing unidentified lines 


TGRS line 
energy, keV 

Count rate 
s' 1 

Transition ID 
(levels, keV) b 

Nominal 
energy, keV 


Comment 3 

63.7 

0.35 

? 65m Ni(63-g.s.) 

63.4 . 

s 


75.6 

0.24 

? 61 Ni(67-g.s.) + K 

75.3 

a 

From 61 Cu(EC). 



? 66 Ga(66-g.s.) + K 

75.5 

a 

From 66 Ge(EC). 

100.2 

0.45 

67 Zn(93-g.s.) + K 

103.0 

a 

From 67 Ga(EC) . 

109.6 

0.22 

19 F(110-g.s.) 

109.9 

s 


246.5 

0.225 

214 Bi(295-53) 

242.0 

r 

238 U series. 



69 Ga(1107-872) + K 

245.2 

a 

From 69 Ge(EC). 

264.9 

0.0031 

? 75 As(265-g.s.) 

264.7 

s 

Weak 

366.1 

0.0055 

? 73 Ge(364-67-13-g.s.) 

364.3 

n 

Three-step transition. 



? 200 Hg(368-g.s.) 

367.9 

a 

From 200 T1(EC). 

384.3 

0.0088 

62 Cu(426-41) 

385.3 

s 


412.8 

0.0082 

228 Th(1432-1023), 

409.5 

r 

232 Th series. 



54 Fe(2949-2538) 

411.4 

a 

From 54 Co (/?+). 



? 198 Hg(412-g.s.) 

411.8 

a 

From 198 T1(EC). 

557.4 

0.0056 

? 62 Cu(548-g.s.) + K 

557.3 

a 

From 62 Zn(EC). 

569.5 

0.0183 

2O7 Pb(570-g.s.) 

569.6 

s 


617.2 

0.0031 

43 Ca(990-373) 

617.5 

a 

From 43 K(/T). 

627.9 

0.0043 

63 Zn(627-g.s.) + L 

628.3 

a 

From 63 Ga(EC). 



? 201m Pb(629-g.s.) 

629.1 

a 

From 201 Bi(EC). 

635.2 

0.0056 

63 Zn(637-g.s.) + L 

638.2 

a 

From 63 Ga(EC). 



63 Zn(627-g.s.) + K 

636.8 

a 

From 63 Ga(EC). 



? n Ge(831-525-198) 

633.1 

n 

Two-step transition. 

645.8 

0.0032 

63 Zn(637-g.s.) + K 

646 , 0 - 

a 

From 63 Ga(EC). 

703.2 

0.082 

? S3m Fe(3040-2339) 

701.1 

a 

From 53 Co(i3+). 



234 U( 1554-852) 

701.6 

r 

238 U series. 

777.9 

0.0056 

228 Th (1168-396) 

772.3 

r 

232 Th series. 

787.2 

0.0073 

? 201 Pb(1415-629) 

786.4 

a 

From 201 Bi(EC). 



234 U(786-g.s.) 

786.3 

r 

238 U series. 

826.2 

0.0036 

? 60 Ni(2159-1333) 

826.1 

a 

From 60 Cu (/?+). 



20 F(823-gs) 

822.7 

s 


851.1 

0.0077 

? 


s 

Energy from simulation. 

947.8 

0.0054 

? 40 Ar single esc. 

949.8 


From 40 K calibration source. 
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Table 4 — Continued 


TGRS line 
energy, keV 

Count rate 
s _1 

Transition ID 
(levels, keV) b 

Nominal 
energy, keV 


Comment 3, 

1005.9 

0.0033 

26 A1 (3074-2069) 

1004.1 

s 




? 48 V(2062-1056) 

1006.3 

s 


1097.7 

0.0044 

228 Th(1153-58) 

1095.7 

r 

232 Th series. 



2O8 Pb(3708-2615) 

1093.9 

r 

232 Th series.: 

1120.6 

0.020 

46 Ti (201 0-889) 

1120.5 

a 

From 46 Sc(EC). 



21 Ne(2867-1746) 

1121 

s 


1163.5 

0.0070 

44 Ca(1157-g.s.) + K 

1161.1 

a 

From 44 Sc(EC). 



62 Ni(2336-1173) 

1163.4 

s 


1190.7 

0.0040 

61 Ni(1185-g.s.) + K 

1193.6 

a 

From 61 Cu(EC). 



234 U( 1237-43) 

1193.8 

r 

238 U series. 

1206.9 

0.00131 

? 200 Hg(1574-368) 

1205.7 

a 

From 200 T1(EC). 

1222.2 

0.0029 

? 24 Na(2563-1345) 

1218.1 

s 


1337.6 

0.0013 

69 Ga(1337-g.s.) + L 

1338.0 

a 

From 69 Ge(EC). 

1378.0 

0.0051 

214 Po(1378-g.s.) 

1377.7 

r 

238 U series. 



57 Co(1378-g.s.) 

1377.6 

a 

From 57 Ni(EC). 

1439.2 

0.0064 

52 Cr(1434-g.s.) + K 

1440.1 

a 

From 52 Mn(EC). 

1609.6 

0.0063 

25 Mg(1612-g.s.) 

1611.7 

s,a 

From 25 Na(,d“). 

1632.1 

0.032 

20 Ne(1634-g.s.) 

1633.6 

a, $ 

From 20 F (/?-). 



23 Na(2076-440) 

1636.0 

s 


1641.2 

0.0015 

? 38 Ar(3810-2167) 

1642.7 

s 


1677.9 

0.0053 

58 Fe(1675-g.s.) + K 

1681.8 

a 

From 58 Co(EC). 

1703.9 

0.0012 

? 2211-2223 keV line 

~ 1700 

s 




complex single escape 




1711.4 

0.0024 

? 73 As(1796-84) 

1712.1 

s 


1893.3 

0.0014 

? 63 Cu (2858-962) 

1895.6 

s 




? 2O3 Pb(2713-820) 

1893.0 

a 

From 203 Bi(EC). 

2092.2 

0.0014 

? 

2092.5 

a 

Line energy from simulation. 

2103.0 

0.00164 

? 2614 keV single esc. 

2103.5 

r,s 


2298.4 

0.0038 

15 N (7566-5271) 

2296.8 

s 




n B(6743-4445) 

2297.8 

s 




? 

2300.5 

a 

Line energy from simulation. 

2373.8 

0.00184 

64 Zn(2374-g.s.) + L 

2375.5 

a 

From 64 Ga(/3+). 
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Table 4 — Continued 


TGRS line 

Count rate 

Transition ID 

Nominal 

Comment a 

energy, keV 

s- 1 

(levels, keV) b 

energy, keV 


2792.2 

0.00076 

? 21 Ne(2794-g.s.) 

2794 

s Weak. 

3393.3 

0.00058 

66 Zn(3381-g.s.) + K 

3391.0 

a From 66 Ga(EC). 


a Key to reaction types: a — /3-decay after activation; s — spallation followed by prompt 
de-excitation; n — (n, 7) or ( 71,71') followed by prompt de-excitation; r — natural radioac- 
tivity. 

b K and L are atomic sub-shell binding energies in cases of electron capture (EC). 


Table 5. A comparison of the simulated and the observed count rate in triangular or “sawtooth’ 5 
shaped features that arise from inelastic neutron scattering in the Ge detector. For the two 
limiting cases of a low and a high estimate of the underlying continuum, see text 


Underlying state, 
keV 

Energy range, 
keV 

Low Cont. Estimate 

High Cont. Estimate 

Sim., s 1 

Obs., s 1 

Sim., s 1 

Obs., s^ 1 

596 ( 74 Ge) 

592-686 

0.66 

0.6 

0.37 

0.28 

692 ( 72 Ge) 

687-760 

0.30 

0.35 

0.09 

0.06 

834 ( 72 Ge) 

831-866 

0.105 

0.053 

0.074 

0.013 

1039 ( 70 Ge) 

1031-1083 

0.166 

0.094 

0.087 

0.070 

1204 ( 74 Ge) , 1215 ( 70 Ge) 

1200-1260 

0.153 

0.109 

0.070 

0.061 
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